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1. INTRODUCTION 


1. 1 Objective of This Document 

The Shuttle Approach and Landing Tests (ALT) will provide 
a large quantity of data on the shuttle navigation system 
capabilities and performance. It is not practical to assume 
that all aspects of the navigation system can be tested during 
the flight itself. For example, it is expected that many 
features of the redundancy management system will not be exercised 
during the ALT. The use of effective post-test processing of 
recorded navigation sensor data can greatly increase the informa- 
tion return from the ALT and from all tests that follow by 
providing the capability for evaluating the existing navigation 
system and for exercising and evaluating alternative versions of 
the navigation system. This can be done without altering or 
increasing the complexity of the actual tests. It is the purpose 
of this document to define cost-effective approaches to such 
post-test data analysis. 

One of the design goals chosen for the development of the 
post-test analysis techniques should be to consider approaches 
that are useful both for the ALT and for the mission data analysis 
of the shuttle after it has become operational. This will greatly 
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aid the performance evaluation and design-modification processes 
during the early operational phases of the shuttle program. In 
this manner, tools developed for the ALT data analysis will 
become more cost-effective by way of maximizing their usefulness 
in other phases of the shuttle development. 

The post-test data analysis activities, as perceived in 
this document, are categorized as follows: 

Quick- Look Processing: Quick-look processing (QLP) , as the 

the term implies, is the preliminary, rapid processing of 
recorded data. This includes such activities as data 
listing and plotting, sensor status reporting, failure 
analysis, and navigation mode monitoring throughout the 
flight. This analysis category allows for a rapid 
determination of test performance as it isolates the 
major problem areas and other unusual events. 

Post-Flight Navigation Processing : Post-flight navigation 

processing (PFNP) refers herein to the execution of the 
navigation functions through the use of flight recorded 
data. Post-flight versions of the operational filters, 
navigation functions, and redundancy management functions 
are used to recreate the navigation portion of the test. 

By altering data and/or navigation software formulations, 
evaluation of all aspects of the navigation software can 


1-2 

INTERMETRICS INCORPORATED • 701 CONCORD AVENUE * CAMBRIDGE, MASSACHUSETTS 02138 • (617) 661-1840 



be tested in an operational environment. In addition, 
alternative formulations can be tested. 

Error Isolation Processing : Error isolation processing 

(EIP) refers to the determination of the magnitudes of 
those error sources that contribute to navigation error. 

By the isolation of the error contributors to a component 
level, problem areas can be identified. In addition, the 
determination of component errors can greatly aid in the 
detailed evaluation of the redundancy management formula- 
tion. 

Each of these three areas of post-test data processing is 
described in this report. Included in the descriptions are their 
respective advantages, disadvantages, uses, design, operation, 
and data requirements. 

1. 2 Format of the Report 

In setting forth the post-test data analysis techniques for 
ALT navigation data, the overall post-test processor requirements 
are described first. In that requirements section, operational 
and design requirements are discussed separately, followed by the 
data input requirements and the requirements on the software test. 

In Section 3, the post-test data processing is described in 
three segments based on the natural test sequence: from quick- 
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look analysis, to post-flight navigation processing, to error 
isolation processing. In each of the three segments, the level 
of detail includes a functional requirements description and a 
description of the software requirements. Section 3 concludes 
with a summary, emphasizing the tradeoffs that must remain open 
and subject to analysis until final definition has been achieved 
in the shuttle data processing system and the overall ALT plan. 

A development plan for the implementation of the ALT post- 
test navigation data processing is presented in Section 4. This 
plan is a necessary adjunct of this report so as to show the 
feasibility of having the system as described herein completed 
in time to support the shuttle approach and landing tests. The 
document closes with a summary and conclusions. 
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2, POST-TEST NAVIGATION SYSTEM REQUIREMENTS 


2 • 1 Operational Requirements 

The ALT phase of the Orbiter development will qualify the 
vehicle for the last and most critical phase of its operational 
use , namely/ approach and landing. To accomplish this test, the 
Orbiter will be carried aloft piggy-back on a modified Boeing 747 
for one or more captive flights. Following these captive flights, 
the Orbiter will be carried aloft and released for a free-fall 
duplication of the approach and landing. 

Because of the cost of each flight, in terms of support 
personnel and systems, the ALT will place a premium on rapid and 
efficient attainment of test objectives. All design features 

should be incorporated in the post— test data processing 
to ensure that all possible information about the navigation 
system performance is extracted from each flight, and that no 
flight test or flight test phase needs to be duplicated because 
of a data processing or analysis deficiency, whether in concept 
or in implementation. 

Because of the need to optimize carefully both the status 
of the Orbiter systems after each flight and the flight test plan 
itself, it is not anticipated that fast turnaround between 
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individual flights will be the rule. Nevertheless, because of 
the need to identify immediately the sources of anomalies, 
discrepancies, and failures, it is anticipated that severe 
schedule pressures will be applied to the test data analysis 
operation following each flight. Following anomaly source 
identifications, it will often be necessary to require hardware 
or software modifications under conditions of limited time and 
approaching deadlines. Original vendors will be called in, 
and there will be redesign, simulation, and repeats of qualifica- 
tion tests. All these activities place tremendous pressures on 
early and efficient problem diagnosis. 

In the light of these and other considerations, the next 
subsection discusses the design requirements for the post-ALT 
data analysis. 

2.2 Design Requirements 

2.2.1 General 

The navigation data analysis system must have a clear, 
simple method for handling all parameters that might vary during 
the entire ALT phase. This means that measurement lists, gains, 
biases, units, must all be alterable under proper project control. 
Alteration must include timely specification in the flight test 
plan, loading, verification, and documentation. The software 
system must be able to run in a pre-flight checkout mode with 
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sample or previous test data so as to provide an unambiguous 
"ready" status in support of the flight. 

2.2.2 Front End 

The front end of the system should be a central file 
maintenance module. As input , this activity should be able to 
accept tapes or other machine-readable forms from as many different 
data sources as possible, with the main design emphasis being 
placed upon telemetry and onboard recordings from the Orbiter 
data processing system (DPS). Additionally, the maintenance 
module must be able to accommodate tapes containing downlinked 
data from the telemetry station, if any; tapes from the acceptance 
checkout equipment (ACE) system in the hangar; and tapes from 
the ground-based tracking organization (s ) . 

The primary design aim of the central file maintenance 
module is to place data of various formats into a data base that 
would be accessible by numerous analysis processes. There is a 
strong design requirement to develop those analysis processes into 
a common higher-order language, as a matter of cost efficiency in 
program development. 

The maintenance module should be self-reporting and self- 
documenting in its performance. It should output a post-run 
analysis that will help with data control and analysis in the 
succeeding phases. Specifically, the post-run output should 
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describe the input tape and format, the output tapes and format, 
the anomalies, gaps, etc., discovered, and the clear failures 
discovered on a first- run basis. 

2,2.3 Blanking and Blending 

Another front end-function involves the ability to prepare 
pseudo-f light tapes for the purpose of simulating the presence 
or absence of selected data categories. This capability would 
enable study of alternative test flight scenarios when a major 
sensor item has failed, is restored, or has been modified. The 
inputs to the blanking and blending module would contain the 
preprocessed flight tapes from the maintenance module; simulated 
time histories of sensor data; and miscellaneous instructions. 

The output would be a new, pseudo- flight tape containing modified 
time histories. This new output tape can be written: 

a. in the same standard format as the input tape if 
substitutions only are involved; 

b. in longer standard format if supplemental time 
histories are added; and 

c. in a shorter standard format if there is blanking 
of certain parameters. 

As in the case of the maintenance module, the blanking and 
blending module should contain a self-documenting listing of the 
input or source tapes, and the modifications to be accomplished. 
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2.2.4 Specific Processor Design Requirements 

The post processors must be able to handle data over the 
entire operational range of the shuttle ALT. The processors must 
also be simple to use and modify and must have sufficient 
commonality with the operational software (variable names , algorithm 
formulation, etc.) to allow easy modification of the operational 
software based on post-test results. The following text presents 
specific constraints that must be placed on the post-test software. 

a. The data analysis processors must be able to handle 

data obtained from the following ALT mission characteristics: 

.(1) Maximum shuttle speed - 1,600 fps 

(2) Maximum shuttle altitude - 70,000 ft 

(3) Maximum shuttle load factor - 2 g 

These constraints have been defined to assume adequate post-test 
processor algorithm formulation and precision. 

b. The post-test processors must be capable of 
simultaneously utilizing the data from the following sensors: 

(1) 3 Inertial Measurement Units 

(2) 3 TACAN receivers 

(3) 3 MSBLS receivers 

(4) 3 Air Data Units 

(5) 2 Radar Altimeters 
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(6) Theodolite Ground Tracking Data 

(7) OB and Radar Ground Tracking Data 

The requirement for simultaneous use of this data in the post- 
test processors is based on the desire to use the post— flight 
processor (see Section 3*2) as a method of evaluating redundancy 
management algorithms* 

c. The post-test processors must be capable of processing 
multiple input tapes* They must also be capable of accepting 
flight-recorded or downlink-recorded shuttle data. This would 
permit maximum utilization of available data. 

d. The post-test processors should be constructed in a 
structured top-down manner. This approach provides increased 
program readability and maintainability. A major goal of this 
approach is to minimize the time required to change portions of 
the post-test processors to reflect changes in output desired, 
candidate navigation or redundancy management algorithms to be 
tested, etc. 

e. The post-test processors should also be constructed 
using the HAL language. This requirement is most critical for 
the post-flight processor. As envisioned (see Section 3.2), the 
post-flight processor contains a post-flight version of the 
operational navigation and redundancy management algorithms. The 
use of the operational computer language, HAL, would facilitate 
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verification of the similarity between operational and post- 
flight versions. In addition, any algorithm changes that are 
made and tested post- flight and are to be incorporated into the 
operational version can be made most accurately by using the same 
computer language for the operational and post- flight computer 
programs . 

2. 3 Data Collection Requirements 

Three major types of data, and a number of minor types, 
must be collected for the design as presently envisioned. The 
major types are the preset navigation parameters and constants, 
the shuttle navigation sensor data recorded throughout the flight, 
and the ground tracking trajectory data. This section presents a 
preliminary list of variables for each of the major data types. 

The few minor, types are introduced within the functional 
descriptions found in Section 3. 

a. The preset navigation parameters and constants are time 
invariant variables such as earth rate, conversion factors, pad 
loaded constants, etc. Many of these parameters are known apriori, 
but others must be recorded at the time of the pad load to ensure 
an exact duplication of the on-board system. A list of these 
parameters (which will be used to generate a parameter file tape; 
see Section 3.2) is presented in Table 2-1*. 


*The data lists presented in Tables 2-1 and 2-2 are based on 
variables presented in the document "Approach and Landing Test 
Level C Requirements for Navigation", Sept. 3, 1974. 
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Table 2-1 


Present Navigation Parameters and Constants 


Variable Name 

Data 

Type 

Variable Description 

ALT MLSAZ 
ALTJTLSR 
ALT_MLSEL 
ALT_TACAN 

Scalar 

Scalar 

Scalar 

Scalar 

Altitude above ref. 
ellipsoid of MLS 
azimuth, elevation, and 
range radars and the 
TACAN station, 
respectively 

AZ RADAR_ 
BEARING 

Scalar 

Bearing from true north 
of MLS azimuth radar 
boresight axis 

BARO DATA_GOOD 

Bit 

' 

Bit set at start- of ^ 
measurement processing 
sequence to indicate 
presence of good baro 
altimeter data : 

BIASACC1 (2,3) 

Vec (3) 

Acceleration bias-vectors 
for the 3 states 

BIAS_ALIGN1(2,3) 

Vec (3 ) 

I MU misalignment biases 
for the 3 states 

bias_azmls 

BIAS__BARO 

bias_dme 

bias__elmls 

b i as _mls range ■ 

BIAS_VOR 

Scalar 

Scalar 

Scalar 

Scalar 

Scalar 

Scalar 

Initial values of the 
sensor biases for MLS 
measurements and TACAN 
measurements and BARO 
measurement 

DELR_RAMP 

Scalar 

Distance over which 
velocity weighting of 
radar altimeter is ramped 
from 0 to its maximum 
value 

e 

Scalar 

Constant natural log 
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Table 2-1 (Continued) 


Variable 

Data 

Type 

Variable Description 

ELLIPT 

Scalar 

Earth Ellipticity 
constant 

EL RADAR 

Scalar 

Bearing from true north 

BEARING 


of MLS elevation radar 
boresight axis 

H_FBAR 

Scalar 

Altitude below which baro 
is inhibited. Automatic 



mode 

HJCNBAR 

Scalar 

Altitude at which baro 
measurements are 
initiated. Automatic 



mode 

H RADALT 

Scalar 

Altitude below which 


5 

radar altimeter is 



■ initiated 

J2 

Scalar 

Value of J2 constant in 
earth f s gravity model 

K RES EDIT 

Scalar 

Scale factor on filter 

, 


mean square residual used 
in filter residual edit 



test 

k_und_wgt 

Scalar 

Measurement underfreighting 
- factor. Multiplies 
BT E B which is then 



added to VAR. " 

K_VAR_AZMLS 

Scalar 

Scale factors on 

K_VAR_ELMLS 

Scalar 

respective sensor 
variances used in LAND 

K_VAR_RN LS 

Scalar 

FILTER residual ‘ edit test. 

K_VAR_RADALT 

Scalar 
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Table 2-1 (Continued) 


Variable Name 

Data 

Type 

Variable Description 

LAT_MLSAZ 

LATJ4LSEL 

LATJ4LSR 

Scalar 

Scalar 

Scalar 

Geodetic latitude of 
appropriate MLS sensor 
radar location 

LATJTACAN 

Scalar 

Geodetic latitude of 
TACAN radar location 

LONG_MLSAZ 

LONG_MLSEL 

LONGJ4LSR 

Scalar 

Scalar 

Scalar 

Longitude of appropriate 
MLS radar location 

LON GJT AC AN 

Scalar 

Longitude of TACAN radar 
location 

R_G°_RAMP 

Scalar 

Value of slant range 
between shuttle and 
landing site below which 
radar altimeter data is 
processed in the filter 

R NAVBASE 
BODYFR 

Vec (3) 

Navigation base location 
in body coordinate 
system 

. R RAD ALT 
BODYFR 

Vec (3 ) 

Radar altimeter location 
in bodv coordinate 
system 

R JTACAN__EF 

Vec (3) 

TACAN position vector 

R TACANT BODYF 

Vec (3) 

Shuttle TACAN receiver 
location in body 
coordinate system 

R MLS ANT 
"'BODYFR 

Vec (3) 

Shuttle MLS receiver loca- 
tion in body coordinate 
system 

R RMLS_EF 

Vec (3) 

MLS range radar position 
vector 
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Table 2-1 (Continued) 


Variable Name 

Data 

Type 

Variable Description 

RADALT_B I AS 

Scalar 

Radar altimeter bias 

RGO_AZMLS_ZERO 
RGO ELMLS_ZERO 

Scalar 

Scalar 

Fixed weighting vector 
range scale factors for 
MLS azimuth and elevation 
measurements respectively 

RE 

Scalar 

Earth equatorial radius 

RHOJ3EAJEEVEL 

Scalar 

Sea level atmospheric 
density 

RJRDAR_BEARING 

Scalar 

Bearing from ture north 
of MLS range radar ‘ 
boresight axis 

r _ GRAV 

Scalar 

Earth radius used for 
gravity oblateness term 

SCALE_HGT 

Scalar 

Alatitude scale height 
for exponential model of 
atmospheric density 

VAR UNMOD ACC 
AT 

Scalar 

Variance of unmodeled 
acceleration times time 
Step. Used in calculation 
of process noise matrix 
for velocity error due to 
platform misalignment 

VAR_ACC_QUANT 

Scalar 

Variance of accelerometer 
quantization error. Used 
in calculation of process 
noise matrix 

VAR_DRIFT_DT 

Scalar 

Variance of error in gyro 
drift compensation. Used 
in calculating process 
noise matrix for platform 
misalignment states. 
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Table 2-1 (Continued) 


Variable 

Data 

Type 

Variable Description 

VAR_AZMLS 

Scalar 


VAR DME 

Scalar 

A-priori sensor measure- 


Scalar 

ment error variances 

VAR ELMLS 

(VAR HI and VAR H2 used 

VAR HI 

Scalar 

for BARO measurements). 

VAR_H2 

Scalar 


VAR_RADALT 

Scalar 


VAR_RMLS 

Scalar 


VAR_V0R 

Scalar 


V RHI 

Scalar 

Upper and lower limits 


Scalar 

respectively of relating 

V__RLO 

wind velocity between 
which BARO updating is 




prohibited. 

C ' 

Scalar 

Difference criteria in 

^AZ 

MLS range, azimuth and 


Scalar 

elevation measurement 

b EL 


change check. 

C R 

Scalar 

respectively 


Scalar 

Difference criteria used 

'^TAG 

to determine if state 
vectors should be 



- 

extrapolated (from T_FILT 
to T SENSORS) to compute 



measurement residuals 

r 

Scalar 

Difference criteria to 

T_MLS 


force MLS azimuth or 
elevation measurement 
incorporation even though 
the measurement is within 
£ AZ or £el °f the previous 
measurement 

y 

Scalar 

Universal .gravitational 

H 


constant for the earth 
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Table 2-1 (Continued) 


11 



Variable Name 

Date 

Type 

Variable Description 

IT 

Scalar 

PI 

a _ ASCENT 

Vec (12) 

■ 

. 

Represents the one-sigma 
values of the initial 
filter covariance matrix 
diagonal elements. In 
local vertical coordinate 
frame (u,v,w) ; used in 
ASNT INITIALIZE 

a ACC 

Scalar 

RMS acceleration bias for 


* 

each of the three filter 
states. a ACC same for 
each of the three accele- 
ration biases of each 
state ... 

a_AZMLS 

Scalar 


cr_DME 

Scalar 

RMS values of particular 

a_ELMLS 

Scalar 

sensor measurement bias 

a_RMLS 

Scalar 


c_V0R 

Scalar 

' 

0J3ARO 

Scalar 

Altitude dependent RMS 
value for BARO measure- 



ment 

a BARO ZERO 

Scalar 

RMS BARO altimeter error 



at zero altitude 

OjJNDJfCT 

Scalar 

Measurement underweighting 
criterion • I 
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Table 2-1 {Continued) 


Variable Name 

Data 

Type 

Variable Description 

T ACC 

Scalar 

Acceleration measurement 

T AZMLS 

Scalar 

MLS azimuth measurement 

T_BARO 

Scalar 

Baro measurement 

T DME 

Scalar 

TACAN DME measurement 

X_ELMLS 

Scalar 

MLS elevation measurement 

t_rmls 

Scalar 

MLS range measurement 

T_VOR 

Scalar 

.. TACAN VOR measurement 

03 ALT_RWTD 
W_ALTD0T 

Scalar 

Prestored gain for radar 
altimeter measurements in 
the altitude and altitude 
rate channel (max value 
of altitude rate gain) . 
Landing phase 

W ALTDOTJWTD 

Scalar • 

Value of altitude rate 
filter gain for radar 
"" -altimeter measurements in 
landing phase 

wE 

Scalar 

Earth rotation rate 

W LAND_AZMLS 
0] LAND_ELMLS 
<i>_LAND_RMLS 

Vec (6) 
Vec (6) 
Vec (6) 

Prestored gain vectors for 
MLS azimuth, elevation, 
and range measurements 
respectively in landing 
phase filter. (Position 
' and velocity updating) in 

the .respective radar 
coordinate systems 
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b. The shuttle navigation sensor data includes the navigation 
parameters that must be recorded throughout the flight. This 
includes sensor outputs with time tags , mode status switches , data 
good flags, and selected computed navigation parameters like the 
shuttle position vector from the navigation state selection filter. 
A preliminary list of those parameters required for post-flight 
processing is presented in Table 2-2. The assumption is made 

that these parameters are recorded at a rate of AT_AVG_G_ENTRY . 

c. The ground tracking trajectory data is an independent 
trajectory as determined from C— band and/or cinetheodolite data 
with associated time reference. It is assumed that that data will 

be presented in a processed form. 

The lists of data collection requirements presented in 
this section define a baseline for the presently-envisioned post- 
test processors. Variations in this list will occur as the 
processor design develops further. In the meantime, this list 
will be useful as a preliminary definition of the data recording 
requirements. 
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Table 2-2 


Shuttle Navigation Sensor Data 


Variable Name 

Data 

Type 

Variable Description . 

AZMLS_EDIT 

Bit 

On-indicates MLS azimuth 
measurement has failed 
filter residual edit test 

AZMLS EDIT 
OVERRIDE 

Bit 

Allows manual override to 
force' MLS azimuth measure- 
ment incorporation even 
though filter residual 
edit test was not passed 

AZMLS_READ 

Scalar 

, MLS azimuth sensor reading 

BARO MODEL OF 
THEJ40NTH 

Bit 

If on / current nav panel 
BARO model is used in 
the filter 



. Data parity bit . . . 



INS attitude from INS 



. . Airdata pressure 



. . . . Temt> measurement . . 

- 


Clock time for sync with 
ground track 

BAR0_EDIT 

Bit 

On-indicates BARO measure- 
ment has failed filter 
residual edit test 

BARO EDIT 
OVERRIDE 

Bit 

Allows manual override of 
filter residual edit to 
force incorporation of 
BARO data 

BARO READ 

. Scalar 

. . . BARO sensor reading 
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Table 2-2 (Continued) 


Variable Name 

Dat a 
Type 

Variable , Description 

DG__AZMLS 

Bit 


DG BARO 

Bit 

Data good discretes from 


Bit 

sensor reads 

DG_DME 


DG ELMLS 

Bit 

• 

DG__RADALT 

Bit 


DG_RxMLS 

Bit 


DG_V0R 

Bit 


DME EDIT 

Bit 

On-indicates DME measure- 



ment failed filter 



" "residual edit test 

DME EDIT 

Bit 

Switch to force TACAN 

OVERRIDE 


range measurement 

* " 


incorporation independent 



of filter residual edit 

— — — - ■ 


test 

DME_READ 

Scalar 

TACAN. DME sensor reading . 

DO_BARO__NAV 

'Bit 

Switch which indicates 

DO MLS NAV 

Bit 

which sensor measurement 

_ — - 


is to be processed 

DOJTACAN_NAV 

Bit 

■ 

DO_RADALT_NAV 



EDIT FLAG 

Bit 

Local variable. On- 


• 

indicates that the sensor 



measurement being 



processed has failed 



filter residual edit test 

ELMLS EDIT 

Bit: 

On-indicates MLS elevation 



measurement has failed 
filter residual edit 

ELMLS EDIT 

Bit 

Allows manual override of 

OVERRIDE 


filter residual edit test 



to force incorporation oi 



MLS elevation measurement 


2-17 


INTERMETRICS INCORPORATED « 701 CONCORD AVENUE • CAMBRIDGE, MASSACHUSETTS 02138 • (617) 661-1840 















Table 2-2 (Continued) 


Variable Name 

Data 

Type 

Variable Description 

ELMLS_READ 

Scalar 

MLS elevation sensor 
reading 

IMU_FAIL 

Boolian 

From IMU RM module. 
Indicates which IMU or 
combination of IMU’s has 
failed 

MANUAL EDIT 
OVERRIDE 

Bit 

Local variable in ENTRY_ 
FILTER allows manual 
override of filter 
residual edit test 
allowing measurement 
incorporation 

-MANUAL BARO 
PRESENT MODE 
MANUAL TAC MLS 
PRESENT MODE 
MANUAL RAD ALT 
"“PRESENT_MODE ' 

Bit 

Bit 

Bit 

Bit 

Allows manual presenta- 
tion of navigation sensor 
data for measurement 
processing 

MANUAL_BAR0 
MAN U AL__T AC_MLS 

manual radalt 

Bit 

Boolian 

Bit 

Switch settings on to 
select BARO, TACAN, or 
MLS , or RADAR altimeter 
respectively 

MLSJTIME_READ 

Scalar 

Time associated with MLS 
sensor reads 

MLS_ID 

Scalar 

Used in identifying MLS 
stations so that pertinent 
information (locations of 
range, az, el, radars, 
etc.) can be obtained from 
stored tables 

RADA LT_E D I T 

Bit 

On-indicates radar 
altimeter measurement has 
failed filter residual 
edit 


2-18 

INTERMETRICS INCORPORATED * 701 CONCORD AVENUE * CAMBRIDGE, MASSACHUSETTS 02138 * (617) 66M840 














Table 2-2 (Continued) 


Variable Name 

Data 

Type 

Variable Description 

RAD ALT EDIT 
OVERRIDE 

Bit 

Allov/s manual override of 
filter residual edit for 
radar altimeter measure- 
ments 

RADALT_READ 

Scalar 

Radar altimeter sensor 
read 

RMLS_EDIT 

Bit 

On-indi cates MLS range 
measurement has failed 
residual edit 

RMLS EDIT 
OVERRIDE 

Bit 

Allows manual override of 
* filter residual edit test 
for MLS range measurements 

RMLS_READ 

Scalar 

MLS range sensor read 

r MID 

Vec(3) 

Shuttle position vector 
from navigation state 
selection filter 

r ONE (TWO , 

THREE) 

Vec (3) 

Current shuttle position 
vectors (3 states) 

SENSJTIME_READ 

Scalar 

MTU sensor read time 

TACAN_ID 

Scalar 

Used to identify TACAN 
stations whose measure- 
ments are to be processed 

T BARO 
T TACAN 
T RADALT 

Scalar 

Scalar 

Scalar 

Time of BARO, TACAN (VOR, 
DME) and 'radar altimeter 
sensor read respectively 

T_COV 

Scalar 

Time tag of filter 
covariance (also used as 
dummy local variable) 
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Table 2-2 (Continued) 


Variable Name 

Data 

Type 

Variable Description 

T_E F _AZMLS 

T _ EF _ ELMLS 

T_EF_RMLS 

Scalar 
‘ Scalar 
Scalar 

Set at start of measure- 
ment processing sequence 
so that a consistent set 
of data is used in the 
measurement processing 
cycle 

T_FILT 

Scalar 

Time associated with 
r FILTl f v^FILTl, etc. 

T_IMU 

Scalar 

Time of IMU read 

T_SENSORS 

Scalar 

| 

Set at start of each 
measurement processing 
sequence; sensor reads 
time 

T STATE 

Scalar 

Time associated with 
current shuttle state 
vectors, r ONE (TOO, THREE 5 
and v _ ONE (TWO, THREE) 

V0R_EDIT 

Bit 

ON-indicated that VOR 
measurement has failed 
filter residual edit 

VOR EDIT 
OVERRIDE 

Bit 

% 

Allows manual override of 
filter residual edit to 
force incorporation of 
VOR measurement data 

V IMU1 CURRENT 
flMU2 , IMU3) 

Vec (3) 

Currently read (T IMU) 
accumulated velocTty 
counts from each IMU 

V HID 

Vec (3) 

Shuttle velocity vector 
from navigation state 
selection filter. Used 
for guidance calculation 
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Table 2-2 (Continued) 


Variable Name 

Data 

Type 

Variable Description 

v ONE (TWO, 

Vec (3 ) 

Current shuttle velocity 

THREE) 


vectors (3 states) 

AT AZMLS 

Scalar 

Time intervals between 

— 


measurements used in time 

AT_ELMLS 

Scalar 

scaling the weighting 

AT RMLS 

Scalar 

vectors for the appropriate 



sensors in the landing 

AT_RAD 

Scalar 

phase filter 

AV IMUl (IMU2 , 

Vec (3) 

Sensed Av from each IMU 

IMU3) 


over AT IMU interval 

AT IMU 




1 



2-21 


INTERMETRICS INCORPORATED * 701 CONCORD AVENUE • CAMBRIDGE, MASSACHUSETTS 02138 * (617) 661-1840 





2,4 Software Checkout 


Development and verification of post-flight processing 
activity requires the use of proper management organization, 
coordination and control techniques. Quality control measures should 
be enforced throughout the development effort. These controls deal 
specifically with test level control, data base control, test case 
generation, design change control,, and software/hardware inter- 
face control. 

The verification of computer programs is greatly enhanced by 
the readability characteristics of the HAL programming language 
and by the use of top-down programming techniques = This allows 
program testing to be conducted in a purposeful, methodical way, 
utilizing a three-level test plan: 

a. Unit Design Verification (Level 0) ~ The objective of a 
Level 0 verification is to validate the lowest-level functional 
units. Each unit is required to compile successfully and is 
further validated by a desk-check of the coding. Input/output 
tests should be utilized, as appropriate, to verify proper control 
flow, numerical results, stability, convergence, scaling, and 
range, etc. 

b. Module Design Verification (Level 1) - Level 1 tests are 
designed and executed as appropriate. The objective of Level 1 
testing is to test related groups or strings of software units 

that represent a higher level of complexity. Applications functions 
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are integrated.. These tests verify inter unit communi cations and 
control, name scope, and parameter passage. Preplanned test cases 
would be used to provide a "calibration standard" aginast which 
results can be evaluated. 

c« System Design Verification (Level 2) - Testing at Level 2 
involves the processing of simulated or real flight data by the 
complete system. This level requires an almost totally-integrated 
system, interfaced with its environment. Such testing reveals 
errors due to conflicts introduced by data-sharing convention 
violations, range of input values, sequencing requirements and 
interphase communications and control. Generated data should be 
evaluated and verified through use of post-processing editing 
software or by visual inspection of the results. 

For the error isolation processor (see Section 3.3) , it is 
required that the simulated nav sensor error sources be contained 
within the estimated + la band about the values of those error 
sources estimated by the EIP. Depending on the trajectory, however, 
some of the error sources will not be observable. This will be 
reflected in large values for the estimated variances of those 
states, thus negating the value of the above suggested test 
criterion. It is therefore important that one sufficient set of 
simulated trajectories, which would make all error sources 
observable at one time or another, be processed in verifying the 
EXP. 
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This three-level test sequence should be followed by acceptance 
tests conducted at the intended residence program facility. These 
tests would validate the interfacing of the system with its 
envisioned environment, allowing the user to confirm compliance of 
the delivered program with its specifications. 
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3. POST-TEST NAVIGATION PROCESS REQUIREMENTS 


The post processing of recorded navigation data can greatly 
improve the efficiency and information return from the shuttle 
approach and landing tests (ALT) . To maximize this return, early 
consideration must be given to the types of processing desired so 
that the processing concepts might have an appropriate impact on 
the data collection processes for the tests. Toward that end, 
this chapter describes post-test processing to a level of detail 
sufficient to outline its impact on the overall shuttle program. 

Three categories of post-test data processing are presented 
in this chapter: quick-look processing (QLP) , post-flight naviga- 

tion processing (PFNP) , and error isolation processing (EIP) . For 
each category, the functional requirements and the software 
requirements are provided. 

Figure 3-1 shows an overview of the interrelationship among 
the three categories of test described in this document. 

3 . 1 Quick-Look Processor (QLP) 

The quick-look processor (QLP) serves two broad post-flight 
navigation data processing functions, each to be carried out 
immediately following the flight test: 1) the QLP serves as the 

central data distributor and in a maintenance function? and 2) it 
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parameters 


Figure 3-1 Post-Test Navigation Data Analysis Process . 




will create and distribute all first-cut navigation evaluation 
reports (i.e., the so-called 24-hour reports). A description 
of the activities within these two functions is presented below. 

3.1.1 Functional Requirements 

The inputs to the QLP are all flight navigation recordings, 
including: in-flight navigation data recordings, both on— board 

and telemetry; preset navigation parameters; and ground based 
tracking data. (See Section 2.3 for the detailed listing of those 
inputs.) The outputs from the QLP are: 1) properly formatted 

versions of the flight data, which will be distributed to naviga- 
tion analysis and design activities, including the two described 
in Sections 3.2 and 3.3 herein; and 2) a series of flight naviga- 
tion summary reports, which will be distributed to any organiza- 
tion having technical responsibility for the performance of 
flight navigation functions.' These two output requirements place 
four functional requirements on the QLP, described in succeeding 
paragraphs. 

3. 1.1.1 Creation and Maintenance of a Central Navigation Data 
File - A central data maintenance function is essential to the 
orderly distribution of complete and readable flight data. Such 
a function alleviates the following potential time-consuming and 
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costly problems for designers and analysts: 

a. inability to locate data, 

b. no time reference, 

c. inability to synchronize two data streams, 

d. improper scaling, 

e. noisy data, 

f . missing data points, 

and a host of other difficulties. In addition, the central 
processor controls the amount of post— f light data processing to 
include only that which is necessary. It eliminates any require- 
ment for data filing, or any other controlled handling, on the 
part of analysts. (This can be accomplished through the maintenance 
of a master copy of all flight data at the central facility.) It 
assures a common time frame for all data streams. It standardizes 
all formats. It helps to assure no duplication of effort. It 
avoids misplacing of data. Also, it minimizes the number of 
constraints on the flight data-recording equipment specifications, 
in that the burden of a common format is on the central facility, 
not on the equipment designers. 

3.1.1. 2 Creation and Distribution of Flight Navigation Summary 
Reports - The primary post-flight operational requirement is to 
enable the many design and analytic talents to work on necessary 
modifications to the navigation system as soon after the flight 
te st as possible . This requires the immediate creation of summary 
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navigation reports to tag anomalies and summarize the performance 
of all navigation functions. And, because of the various and 
different perspectives of the groups concerned with navigation, 
there is a requirement for a number of different summary views 
of the flight (e.g., sensor people will be concerned with readout 
performance, computer people with mode switching, and so forth). 

The second QLP requirements is, therefore, a set of 
summary reports, each with a different perspective of the flight. 
That set of reports represents an organized statement of the 

j 

flight navigation performance. 

These reports should be created within 24 hours of the 
flight, or at least within 24 hours of the availability of 
recordings, in order to expedite analysis of the flight. 

3.1. 1.3 Preliminary Analysis of Unusual Events - The above- 
mentioned summary reports will, among other things, flag unusual 
(unexpected or undesired) events like failures. As a third 
requirement on the QLP, the post-flight data processing specifica- 
tion should include a list of supplementary processing activities, 
each activity within the list dictated by the nature of a 
particular unusual event. For example, a particular sensor failure 
might dictate a processing requirement on that sensor's redundant 
equivalents, or on the ground tracking data, or on its calibration 
data, and so forth. As distinguished from the fixed summary 
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reporting above, this QLP function exercises a set of causal 
requirements (that is, given an unusual event A, accomplish 
activity B) . For this part of the QLP, only those unusual events 
that require immediate action should be included. 

3. 1.1.4 Creation of Modified Quick Looks - The outputs from the 
above three QLP activities will be, for every flight, to a 
variable extent, less than complete. There will always be unfore- 
seen events which will not be covered by the post-flight 
processing specification. Nevertheless, the post-flight data 
processing specification must direct at least the mechanics for 
expediting action on such events. (The specification cannot, 
of course, dictate the particulars of action against unforeseen 
•events.) For the purposes of this document, we assume the 
existence of such a module as the fourth requirement on the QLP. 

The software requirements section that follows contains little 
in the way of particulars, which are left as an exercise for the 
specification. 

Note that there will be a learning curve applied to the 
QLP: action taken in the modified quick-look area on any flight 

should be covered as a specification revision for unusual event 
analyses on the next flight. 

3.1.2 Software Requirements 

The QLP routines required to accomplish the functional 
requirements presented in Section 3.1.1 are introduced in this 
subsection. 
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In Figure 3-2, an overview of the QLP shows it to be 
composed of six routines. For our purposes, herein, "routine” 
is not restricted to mean one distinct physical program. In 
fact, in most instances it is more than one. Although the function 
of any routine might be fixed, there is often a requirement to 
develop a number of programs to accomplish that function, either 
because the software maintenance requires a modular construction 
(because of ever-changing requirements) , or because the physical 
characteristics of the various inputs (e.g., tape, disc, blocking, 
format) dictate several mechanical versions. 

Central to the QLP is the navigation data master file 
(see Figure 3-2). To create that file, two routines are 
required: the Duplicate Flight Recordings (DFR) , and the Purge, 

Format, and Catalog Routine (PFCR) . The DFR simply creates a 
version of the flight test data that can be dedicated to naviga- 
tion analysis alone. This assumed that navigation and non- naviga- 
tion data will be found on the same flight recordings. It also 
assumes that the tape duplications are not accomplished by other 
agencies. The second routine, the PFCR, is that multi-purpose 
activity (probably consisting of many programs) that standardizes, 
formats, catalogs, sorts, and tabulates the flight data. 

From the master file, four routines are required to create 
QLP outputs. The Duplicate Navigation Data (DND) activity acts as 
the central data distributor. Analysts of flight data will be 
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Figure 3-2 Routines Comprising the Quick-Look Processor 









given data dedicated to their use alone. Therefore, they should 
be presented duplicates of the master file copy. The Create 
Flight Navigation Logs activity creates a set of summary reports, 
in detail sufficient for first-look analysis. The Event Analysis 
Routine (EAR) and Modified Quick-Look Routine (MQLR) create 
special analyses and reports for unusual events; the EAR for events 
covered by the specification, and .the MQLR for those not covered 
by the specification. 

In the following paragraphs, the software requirements for 
each of these six routines are covered. 

3. 1.2.1 Duplicate Flight Recordings ( DFR) - This activity may 
or may not exist within the QLP , depending upon the inter- 
relationship between the navigation post processing and non-navi- 
gation post processing activities. Assuming that navigation and 
non-navigation data will exist on one or more of the flight 
recordings, there will be a requirement for duplicating those 
recordings (or partitioning them, which would probably be a 
solution technically less practical than duplication) so as to 
have a dedicated version of flight data for navigation post- 
processing needs. A high level of planning might judiciously 
place that duplication function elsewhere. Until that is specified, 
we herein assume that the function belongs to the QLP. The soft- 
ware requirement is to have the capability to create "exact" 
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duplicates of the flight recordings. For some, there may be a 
temptation to accomplish additional processing during this 
duplication process. This should not be done. Experience m 
most data processing activities shows that an original version 
of data should always be retained, because after data is 
processed, it is difficult to retain a picture of which portion 
of the processed output is raw and which has been altered. 

With an original on file, one can always recover that distinction 
when necessary. Nevertheless, any duplication process is 
difficult to accomplish exactly ; all duplicating processes have 
their own noise. In any event, care should be taken to 
accomplish the duplication as nearly as possible. 

3. 1.2.2 Pu rge, Format, and Catalog Rou tine (PFCR) This 
routine constitutes the most involved set of programs within the 
QLP, due to the large number of functions to be accomplished and 
to the lack of structural relationships among the functions. The 
varied characteristics of the inputs necessitate a "shopping list' 
of standardizing and cataloging processes to enable the inter- 
changeability and cross-referencing of data necessary for the 
reporting processes. The list of software processing require- 
inents is as follows: 

Format Standardization: The input data is of various mechanical 

and logical input formats. The need to cross-reference data 
from the various sources dictates a standardization of format for 
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the master file. This creates a range of process requirements 
including: 


• D a ta Storage Medium Standardization - The input data 
probably will be physically contained in several 
different machine- readable forms. The PFCR should 
translate the data to a standard form, probably tape 
or disc, as dictated by the available processing 
equipment. 

. Data Packaging - The master file is made up of n 
packages of data, with the word "package" used 
here to mean a large collection of data of a 
particular type (e.g., a stream of accelerometer 
outputs). The packaging can be accomplished logically, 

. or physically? the choice would be dictated by 
tradeoffs directed toward efficient retrieval practices. 

. Data Blocking - It is assumed that the processing equip- 
ment and practices dictate a standard quantum (i.e., 
"block") of data. A package would be partitioned into 
n such blocks within the PFCR. 

. Data Logical Record Layout - All data should be trans- 
formed to a common "record" layout, such that a block of 
data equals n records. The specifics of record layout 
would be dictated by efficient retrieval. 
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. Data Addressing - For a tape-oriented system, data 


addressing is a trivial problem. For a disc-oriented 
system, addressing can be used for efficient retrieval. 
However, addressing logic development can be carried 
to extremes. The amount of versatility required for 
accessing through sophisticated addressing techniques 
must be dictated by the output requirements. Care 
should be taken not to exceed those requirements. 

. Data Labeling - Data groups should be properly labeled? 
that is, all packages, all blocks, and in selected 
cases, some records. 

Data Synchronization and Calibration : It is assumed that the 

input data will be of various units and will be referenced to 
several different clocks. A need for conversion to common frames 
of reference creates the following software requirements: 

. Time Conversions to a Master - The different .time 
references for each recording should be brought to a 
common reference if precision requirements so dictate. 
This can be accomplished by event commonalities 
between recordings, like mode switching, failures, etc., 
or by pre- flight clock calibrations. . 
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• Calibration to Similar Units - The instrument data 
strings will have different scaling and biases (and 
perhaps other sensitivity variances, such as 
temperature) . The PFCR should convert instrumentation 
data from various sources to the same units, to that 
extent dictated by the analysis requirements. 

. Calibration to Similar Dimensions - Data from the 
different sources will be attempting identical state 
measurements (e.g., velocity, position) but will be 
of different dimensions, differing only by time 
derivatives or integrals. The PFCR may be the practical 
place to accomplish the integration or differentiation 
to achieve a common dimension. The alternative process 
location would be the analysis routines that require a 
common dimension. A study of the analysis specifications 
will dictate the choice. 

. Line-up of Data Streams - It may be practical to line up 
data streams that fall into a particular class, e.g., 
all gyros, all velocity measurers, etc. The PFCR would 
be the practical place to accomplish this. 

. Bit Stream Conversion - When data is in the form of bit 
streams, it may be advantageous to convert those bit 
streams to strings of digital bit counts. 
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Data Filtering : Data streams will often contain undesirable 

or unnecessary components. The PFCR should be assigned the task 
of removing those anomalies. The following requirements are 
placed on the software: 

• Data Roundoff - Where large packages of data contain 
either unnecessary or meaningless digits, there will 
be an incentive to round-off to the most significant 
data. 

• Data Sample - When it is reasonably certain that the 
data is of much finer grain than necessary, cost savings 
can be achieved in analysis by sampling at less 
frequent intervals, or by summarizing data within 
intervals . 

. Noise Abatement - The removal of non-navigation data 
noise should be accomplished, but only where it is 
reasonably certain that only non-navigation noise is 
being purged. 

Data Catalog : The PCFR would be the logical place to "inventory" 

flight recordings. There are two requirements: 

. Data Completeness - The existence of any loss of data 
during the flight should be recorded and output from 
the PFCR. 

, Data Totals - The statistics of data packages in terms 
of size, granularity, precision, and type should be 
recorded and output. 
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Data Security : Only one trivial software requirement is necessary 

for security purposes, and that is to have more than one copy of 
the PFCR output available, as a protection against fire or other 
damages to the master. 

3.1.2. 3 Duplicate Navigation Data ( DND) - The maintenance of a 
central master file is the best method for sustaining the integrity 
of the total navigation package. Given that the central file is 
maintained as n distinct packages of data (perhaps n separate 
types), the only software requirement for universal distribution 
purposes is a simple duplication process (the DND) . In order to 
have an additional capability for copying any specified subset 

of data imaginable, one would require the DND to have an extract 
logic capability. After thorough consideration is given to the 
organization of the master file, when creating the specification, 
it will probably be found that the versatility of an additional 
extract logic routine will not be required. 

3 . 1.2. 4 Create Flight Navigation Logs ( CFNL) - Flight summary 
reports will be required in order that the various navigation 
design and analysis interests might expedite any required remedial 
design activity for subsequent flights. This remedial activity 
must be begun as quickly as possible. The Create Flight Naviga- 
tion Logs (CFNL) routine (along with the EAR and MQLR) should 
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provide to the analysts, within 24 hours of flight, an overview 
of flight navigation performance. 

The input to the CFNL is all of the data streams within 
the master file; the output is a set of summary navigation reports 
presenting various perspectives of thos data streams. To accomplish 
this output, the following are software requirements; 

Navigation Event Catalog : The flight navigation scenario is given 

by the "navigation state" versus time or versus flight events. 

The "navigation state” is the set of measured and computed para- 
meters that best characterizes the navigation system. At the 
minimum, it will be all instrument recordings (which have been 
formatted, calibrated, and synchronized) and the navigation 
velocity (V) /position (R) state vector. In addition, it will 
include a characteristic set of intermediate navigation computa- 
tions, and any co mputed guidance velocity/position state vectors, 
(i.e., required velocity/position).- 

The assumption at this time is that any specific group of 
analysts will be interested only in a selected subset of the 
navigation state/ which always includes the velocity /position state. 
Also each group would have its own time f rame-of-interest. 

Therefore, the software requirement of the CNFL is to create the 
navigation state vs. time and/or vs. flight events stream, and to 
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partition that package of data according to the type of reader? 
where at least one of the reports would be a summary of the 
whole navigation state vs. time and/or vs. event. 

The events that will dictate navigation state measure- 
ments are: 

. all mode switches , 

. all failures, 

. all Redundancy /Management switches, 

. any unexpected event, and 
. a specified set of time intervals. 

For all types of events there need not be a complete navigation 
state measurement. For example, a Redundancy Management switch 
might require the reporting of only the instruments involved. 

Data Statistica l Calculations: The CFNL should also be required 

to characterize quantitatively each component of the navigation 
state. A minimum set of calculations for such a characterization 
includes: 

, Specification Check - Each critical data stream 
is compared with its specification, and deviations 
are reported. 

. parameters Statistics - The meaningful statistics of 
any parameter are the statistics of the difference 
between that parameter and a set of standard references, 
including 
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* spec limits 
. expected values 

. another parameter that measures the same 
dimension 

. null 

The CFNL should have the capability for creating 
the distribution of any of these relative values. 

. Recording Data Statistics - from PFCR 

. Noise Statistics - any clearly identified noise 
should be reported. 

The statistical calculations will include, additionally, 
a set of specialized data statistics as dictated by analysis 
requirements , and as outlined in the post-flight data processing 
specification. 

3. 1.2. 5 Event Analysis Routine (EAR) - Certain analyses will be 
required only when triggered by unusual events. For example, 
when a particular data stream is suspect, like an accelerometer 
output, this would dictate the comparison of that stream with 
other measurers of the same parameter, like the redundant 
equivalents of the suspect instrument. Because of the complexity 
of the shuttle instrument package, the number of such causal 
analyses can be quite large. The number of types of software 
analyses will, however, be somewhat smaller. The post-flight 
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navigation data preprocessing specification will dictate the 
specific set of such software analysis capabilities. Following 
is the set of requirements that should be expected in that 
specification: 

Data Plotting: When an unusual event is found in any data string, 

the problem is to find the cause or reasons for that event. The 
obvious recourse is to "plot" the string in which the event 
occurred against other strings that might give insight. (To 
"plot", in this case, is to create data pairs that can be 
presented in graph form or in any manner.) The software require- 
ment is to be able to plot any navigation state parameter against 
any other state parameter. Examples of the reference parameters 
are: the velocity /position state; the redundant-equivalent 

string of data; and the ground-track equivalent string. 

Plotting can. also include a finer grain of the analysis 
triggering parameter against time or against events. 

Data Analysis: Any data string or the map of two or more data 

strings can be characterized by certain statistical parameters. 
Examples constituting the minimum software requirement for the EAR 

are : 

. correlation coefficients, 

. means , 

standard deviations, and 
. curve fits (e.g. , least square) . 
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Complementary Data Reporting: There exists within the master 

file a complete set of navigation data for all flights* To aid 
in the diagnosis of unusual events, that file can be used to 
create a "Diagnosis Book" for each class of unusual events* 

Included in that book would be such things as: 

. calibration data from this flight and all previous 
flights 

. instrument repair history, 

. failure history, 

. configuration history, 

and any other data deemed as a necessary input to diagnosis. 

3. 1.2. 6 Modified Quick-Look Routine (MQLR ) - The Modified Quick- 
Look Routine (MQLR) anticipates that thorough foresight at the 
time of specification writing is not possible. There will always 
be a need for quick analyses that are not covered in the preceding 
QLP routines. What the specification can cover, however, is the 
scope of MQLR activity, including, for example, the types. of 
analyses that will be accomplished; level of manpower/machine 
power available; etc. It is expected that the MQLR software 
capabilities will be variations on existing capabilities like 
the EAR routines. It is also expected that MQLR requests will 
expedite revision action against the post-flight navigation data 
processing specification. 
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3.1.3 QLP Summary 

The quick-look processor, as described in the previous 
sections, enables a rapid and orderly distribution of information 
concerning flight navigation performance. Significant cost 
savings will be the result of implementing the QLP as described 
because of: the minimization of waiting on the part of design 

and analysis activities? the elimination of duplicate efforts; 
and the reduction of post- flight testing to include only that 
which is absolutely essential. 
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3.2 Post -Flight Navigation Processor (PFNP) 

The post- flight navigation processor provides the capability 
of executing the redundancy management and navigation functions 
subsequent to flight. A general description of this processor and 
its uses is contained in the following paragraphs. 

3.2.1 Functional Requirements 

Three types of inputs are needed for the post-flight 
navigation processor. These include: preset navigation parameters 

and constants , time history of shuttle navigation sensor data 
recorded throughout the flight, and ground tracking trajectory 
data. Specifics as to the data types are given in Section 2.3. 

The input data is processed by a post-flight version of 
the operational redundancy management (RM) and navigation functions. 
The baseline version of the PFNP should contain an exact duplica- 
tion of the on-board RM and navigation software (filters, sensor 
selection logic, navigation mechanization, etc.). The baseline 
PFNP should be able to reproduce the navigation portion of the 
actual flight using the recorded data. 

The output of the PFNP is a time history of the position, 
velocity, and attitude of the shuttle and filter-determined 
variances of the errors in these quantities. In addition, measure- 
ment statistics for each sensor can be compiled. The trajectory 
output of the PFNP is compared with the trajectory prepared by 
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ground tracking equipment to achieve a measure of the performance 
of the operational software. 

Naturally, the post- flight reproduction of the actual 
test by the baseline processor does not in itself seem too important. 
If the operational software functions properly, such data would 
be available from the test. The real value of the PFNP lies in 
the fact that many alternatives can be exercised using real data 
collected in an operational environment without repeating the 
actual test. Hence, the major uses for the PFNP are: 

Navigation Algorithm Improvement Studies - Changes in the filters 
or other navigation algorithms can be made, and the impact of these 
changes can be tested in the operational environment. Comparison 
of the altered PFNP output with the trajectory determined by 
ground tracking equipment can be used as a performance index. 

Redundancy Management Testing - It is not expected that the ALT 
will exercise all aspects of the redundancy management function. 

By altering the test data, these capabilities can be tested in an 
operational environment. Also, alternative RM algorithms can be 
tested. This can include changes to the failure detection and 
identification (FDI) algorithms. 

Subsystem and Software Error Detection - Post-flight processing 
of the form described in this section is useful for the identifica- 
tion of subsystem and software errors. Three ways in which such 
errors may be identified are: by observing systematic deviation 

of PFNP and ground tracking trajectory solutions, by using 
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expanded state filters, and by comparing sensor measurements 
not used in the PFNP solution with this solution. The latter 
is performed by executing the PFNP using all but one set of 
sensor data. This use of a post-flight processor has been shown 
to be both practical and successful with the CIR1S system*. 

The advantages of developing a post-flight processor for 
the shuttle ALT test are obvious from its potential uses as 
described above. The disadvantage might be the development and 
operational cost associated with this technique. However, this 
cost is minimal with respect to the cost of the test program 
or even one test. 

Compared to pure simulation, post-flight processing as 
described in this section has the advantage of using operationally- 
obtained real data. Modification of a simulator to perform 
this function is often not possible since the post-flight 
processor timing must be governed by the recorded data timing 
instead of by the simulator software. Some of these differences 
will become more apparent in the software description section. 

3.2.2 Software Requirements 

The suggested design of the shuttle post-flight navigation 
processor consists of three separately-operable software modules. 
T hese modules are required to perform the following functions: 


*Radi.o transponder survey errors, doppler radar anomalous performance, 
and INS mechanization equation errors were all found using a CIRIS 
post-flight processor. 
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a. flight data tape editor (FDTE) 

b. parameter file tape generator (PFTG) 

c. navigation function processor (NFP) 

(1) NFP control 

(2) job control card processing 

(3) redundancy management and navigation processing 

(4) ground tracking comparison module 

The system configuration is shown in Figure 3-3. The following 
subparagraphs specify the required operational functions. It 
should be noted that much of the text that follows is a function 
of the method of generating the formatted navigation data. In 
general, multiple tape input will be considered. 
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3. 2. 2.1 Flight Data Tape Editor (FDTE) - This preprocessor j s 

module processes the formatted navigation data tapes and creates 

a new edited flight data tape. This consists of processing the 
flight data to a state of completeness by simulating failed 
sensors and deleting data anomalies. The flight data tape editor 
module has been designed as a separate program for several reasons. 

In many situations, many different PFNP runs can be made while 
the data editing needs to be done only once. Also, in general, 
the same edited data tape can be used for both the PFNP and error 
isolation processor (see Section 3.3). In addition, the flight 
data tape editor provides a standard interface for use in 
generating test tapes used in program checkout. 

3.2.2. 1.1 Source and Type of Inputs - The inputs to the flight 
data tape editor are a subset of the formatted navigation data 
tapes described in Section 3.1.3. Those tapes will have been 
specifically designed for the purposes of the post-test data 
processor immediately following the flight. The master of these 
tapes will be kept in a central file, eliminating a need for 
rigorous control during post-flight navigation processing, 

3. 2.2. 1.2 Destination and Type of Outputs - The output of the 
FDTE should consist of error messages and an edited flight data 
tape or tapes (EFDT) . Output data is generated once per execution. 
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Listable output of messages and processed data should be included 
as an option. 


3.2.2. 1.3 Information Processing - A simple conceptual flow 
diagram of the flight data tape editor is presented in Figure 3-4. 
The job control functions cards are the first processed to 
produce a tape identifier. Error messages for the function cards 
should be included in this processing. Processing is terminated 
if function card errors are found. Next, the navigation data 
is read a record at a time. Errors in the records are identified 
and listable output messages provided. Such errors might include 
parity errors, sync errors, record length errors, and word or 
bit pattern errors. If no errors are found, the record is 
processed. The edited data is then output to a tape and listed, 
if desired. The next record is then read. At the last tape 
record, a check for multiple tapes is made, and, if more than one 
tape is available, a new tape unit is assigned, and processing 
continues . 
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3.2. 2.2 Parameter File Tape Generator - This preprocessor 
module creates the fixed data base for the post-flight navigation 
processor. The data base contains navigation parameter and 
constants as specified in Section 2.3. The sources of data may 
be punched cards, paper tape, etc. The parameter file tape 
generator (PFTG) has been designed as a separate program since 
a common tape input provides a more convenient method of data 
base input to the PFNP. As will be seen in Section 3. 2. 2. 3, the 
inclusion of data base alteration options in the navigation 
function processor (NFP) allows simple data base changes without 
the requirement for re-execution of the PFTG. 

3. 2.2.2. 1 Source and Type of Inputs - The input to the parameter 
file tape generator (PFTG) is of the following three record types: 

. job control functions 
. data 

. job termination 

The job control functions perform housekeeping tasks and 
ultimately result in a tape identity record. This identifier 
should contain such items as tape generation identifier (PFTG) , 
tape identifier, date, originator identification, etc. 

The data records provide for the numeric va3.ues that make 
up the items in the data base (see Table 2-1). The record design 
should provide the name of the value and the value itself. This 
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syntax allows both ease of data input and verification. No data 
record order is required since the envisioned design contains a 
sort routine as part of the baseline. This capability allows 
easy data input and reduces possible errors. 

The job termination records cause the PFTG processor to 
stop reading records and to generate the parameter file tape 
(PFT) . 

3. 2. 2. 2. 2 Destination and Type of Output - The outputs from the 
PFTG module consist of a listing of the contents of the PFT, the 
tape itself, possible errors in the job control functions, data 
variable name errors, and data values which are missing. 

3.2.2. 2.3 Information Processing - A conceptual flow diagram of 
the parameter file tape generator is shown in Figure 3-5. The 
job function records are processed in the same manner as in 
Section 3. 2. 2. 1.3 to produce a tape identifier label. Data records 
are processed by comparing the inputted variable identifier to 

a pre-stored table. If a match is found, the inputted value is 
stored in the appropriate location. During record processing, a 
count is kept of errors encountered in reading the records. 

Presence of any errors will inhibit the generation of the parameter 
file tape; an error summary would be sent to the listable output 
device, and the job would terminate. 
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. If no errors are encountered, a label is written on the 
tape, followed by the prestored table of variable identifiers 
and the table of their corresponding values as modified by the 
data input cards , and the tape is rewound. The contents of the 
parameter file are then listed for record-keeping purposes, and 
the job terminates. 

i 

3. 2. 2. 3 Navigation Function Processor (NFP) - This main module 
of the PFNP recreates the execution of the onboard redundancy 
management and navigation software by operating on the flight 
recorded data. The parameter base is made available through a 
parameter file tape (PFT) ; the contents of this parameter file are 
modifiable under user control. Processing is directed through 
user— generated job control cards. In order to offer maximum use- 
fulness, the NFP provides a large degree of flexibility through 
user options like selection of processing start and stop times, 
selection^ of sensor data to be incorporated by the filter, frequency 
of sensor data processing, etc. 

3.2.2. 3.1 Source and Type of Inputs - Inputs to the NFP consist 
of an edited flight data tape, a parameter file tape, and a job 
control deck. The EFDT is generated by the RFDTE pre-processor 
(see Section 3. 2. 2.1). The PFT is generated by the PFTG pre- 
processor (see Section 3. 2. 2. 2). The job control deck can contain 
cards of the following types: 

, function 
. data alteration 
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. user option generation 
. listable output control 
. termination 


The function card performs administrative duties, leading 
to identification labels to be placed in any tapes and/or listable 
output to be generated during the run. 

The data alteration cards allow modification of the data 
base values stored in the parameter file without requiring its 
regeneration by the PFTG program. 

The user option generation card implements user processing 
directions such as sensor data to be processed, selection output 
tape generation, etc. 

The listable output control card is used to specify the 
frequency and type of data to be listed during the run. 

The termination card indicates the end of the job control 

deck. 

3. 2. 2. 3. 2 Source and Type of Outputs - NFP outputs consist of 

data tapes for further processing or for automatic comparison 

with outside references, and listable output for quick-look analysis. 

All output tapes are labeled with identifiers provided 
through the job control functions, in order to facilitate correla- 
tion of tape data to the conditions of the run in which it was 
generated. Examples of candidate data for tape output include the 
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filter estimate of shuttle position and velocity for comparison 
to cinetheodolite and/or C-band radar reference trajectories, and 
to the onboard computed estimates for verification of compatibility 
of the operational and post-flight software and evaluation of the 
effect of limited computer word precision on the performance of 
the onboard estimator. The Kalman filter state vector and 
covariance matrix allow insight into the filter performance by 
observation of the developed cross-correlations between state 
elements and by comparison of the state variances to the actual 
level of the filter estimation error as generated from ground 
tracking reference trajectories. An important index of filter 
performance is the degree of consistency between the measurement 
residuals (difference between measured and estimated value of the 
measurement) and the filter-computed la error in the estimated 
measurements. Similarly, means and standard deviations can be 
computed from measurement residuals to evaluate the approximation 
of the actual residuals to the zero mean random residual measure- 
ment error assumed by the filter and the appropriateness of the 
selected la error level for the measurement. These processed 
residuals give an indication of the adequacy of the error model 
being used by the filter. 

Listable output would consist of reproduction of the job 
control card deck, listing of the operating data base, error messages 
and diagnostics, and a subset of the tape output data to allow 
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quick-look analysis of the NFP performance before embarking on 
further data reduction. 

3.2.2. 3.3 Information Processing - A simplified flow diagram 
for the navigation function processor is presented in Figure 3-6. 
Following an initialization step, in which administrative functions 
like I/O unit assignments, etc., are accomplished, the parameter 
file ’tape is read. The program then reads the job control cards, 
identifying their types: function cards are processed to set up 
label information for output; data cards are read to modify 
contents of the parameter file data base; generate cards are read 
to set flags which implement user options to control- processing; 
listable output cards are read to determine run pertinent require- 
ments. Reading of a termination card completes processing of the 
job control deck. If errors have been encountered in processing 
the cards, the run terminates following printout of an error 
summary. 
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If so selected, an updated PFT can be generated using 
the data base modifications input through data alteration cards. 

The data base is then listed for record-keeping. Parameter units 
are then converted to a consistent set (°/hr to rads/sec, for 
example) . The program then searches the EFDT for the selected 
processing start time by reading records until the record time 
gxceeds the desired start time. The Kalman filter state vector 
and covariance matrix are then initialized and output to tape 

and the listable output device. 

If an end-of-file has not been reached on the EFDT, and 
the desired stop time has not been exceeded, the program settles 
into its main loop. The next EFDT record is read. If a data 
gap exceeding the filter time step constraint is encountered, an 
interpolated data set is generated for the filter's state and 
covariance matrix update. The data is then processed by a replica 
of the on board- redundancy management and navigation software. 
Output is then made to tape and printer, and the program recycles 
through its main loop. 

On encountering an end— of- file on the EFDT, or on exceeding 
the desired processing stop time, a final output is made, and the 
program terminates. 


3-42 

tNTERMETRiCS INCORPORATED • 701 CONCORD AVENUE * CAMBRIDGE, MASSACHUSETTS 02138 * (617) 661-1840 



3.2.3 PFNP Summary 

The concept and formulation of the post-flight navigation 
processor described in the previous sections provide not only a 
post- flight navigation data reconstruction capability, but also a 
means of altering and testing with flight data the redundancy 
management and navigation functions. These capabilities offer 
a significant means of maximizing test data returns while 
minimizing test operational costs and risk. If formulated using 
the guidelines of Section 2, modification and checkout of the 
operational software based on post-test analysis can also be 
simplified. 

3 . 3 Error Isolation Processor (EIP) 

The error isolation processor (EIP) serves a twofold purpose: 
on the one hand, by estimating the operating error sources, it 
provides a capability to evaluate the performance of the redundancy 
management logic; on the other hand, it provides the capability 
of identifying those error source, or groups of error sources, 
whose contribution to the system error is critical to the 
performance of the shuttle navigation function. A general 
description of this processor and its applications is given below. 

3.3.1 Functional Requirements 

Two input data types are required for the error isolation 
processor: modeling parameters, and time history of shuttle 

navigation sensor data recorded throughout the flight. In addition, 
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ground tracking trajectory data can be used to verify the 
performance of the processor, or as an external reference source 
in the estimation process. 

The input data is processed by a high-order Kalman filter, 
modeling all of the significant navigation sensor error sources. 

In order to make optimal use of all the available information, the 
filter model could include the significant error sources of each 
of the three shuttle IMUs. Given a conservative estimate of at 
least 40 significant error sources for each IMU, the required 
state dimension, with its related computer storage and computa- 
tion time requirements, rapidly becomes unmanageable. Furthermore, 
the errors from the different IMUs can be expected to be highly 
correlated, given that they are driven by the same environment 
(gravity errors) and that they maintain a fixed relative orienta- 
tion. Therefore, their contribution to improved estimation of 
each other's error sources will be small. It therefore seems 
more cost-effective to optimally blend each IMTJ’s outputs with 
the external reference data one at a time. The same processor 
can still be used by providing for user selection of the IMU data 
to be processed. The state dimension is reduced to a manageable 
level ( ~ 5 0 states) with corresponding savings in program run time. 
Note, however, that a source of information has been neglected 
(i.e., the other two IMUs) and that, therefore, the processor 
performance is somewhat degraded versus the theoretically optimal 
achievable estimation accuracy. 
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The output of the error isolation processor is a time 
history of the estimates of vehicle position velocity and 
attitude and of the modeled nav sensor error sources, and the 
filter- computed covariance of the estimation errors. The trajectory 
output can then be compared to the ground tracking equipment 
trajectory to provide confidence in the filter performance. 

The desired output from the processor are the estimates 
of the operating navigation sensor error sources. Two major uses 
for this data are presented below: 

a. Redundancy Management Software Evaluation - Determination 
of the IMU error sources provides the basis for RM evaluation by 
establishing whether a failure had indeed occurred when the soft- 
ware so indicated , or conversely, whether the software failed to 
react in the face of an actual component failure. On this basis, 
the RM design can be validated or appropriate modifications 
justified. 

b . Identification of Error Sources Critical to Navigation 
Performance - An error source, or a linear combination of error 
sources, is observable to the extent that is contributes to the 
measured error. Thus, in the event that the navigation filter 
estimates are not consistent with external measurements, resulting 
in degraded filter performance, the error isolation processor 
provides a ready tool for the identification of error sources, 
modeling errors, etc., causing the poor performance. With its 
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high-order state, and therefore greater error dissemination 
capability, it allows accurate determination of the error source 
characteristics and permits clear-cut testing of competing 
hypothesis . 

This processor also offers the possibility of using 
smoothing techniques to achieve improved estimation accuracy. 

Smoothing (in the sense of Fracier) represents the combination 
at any point in time of optimal filters processing the data forward 
and backward in time. It thus provides, at any point, the optimal 
estimate derived from all the available data, with significant 
improvement of the estimation of time-varying errors (for constant 
errors, the final forward filter estimate represents the result 
of processing all the available data, and is, therefore, the optimal 
estimate) . The benefits to be provided by smoothing depend then 
on the extent to which the shuttle nav sensor error sources are 
expected to be time-varying. In general, most I MU error sources 
are considered to be constants, and smoothing would not then provide 
any improved estimation accuracy. 

3.3.2 Software Requirements 

The suggested design of the shuttle error isolation processor 
consists of three separately operable software modules: 

a. flight data tape editor (FDTE) , 

b. parameter file tape generator (PFTG) , 

c. error analysis processor (EAP ) . 
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The first two modules follow the same design approach 


as those described for the PFNP; in fact, a generalized version 
of these modules could be developed to serve as interfaces for 
both main processor (PFNP and EIP) . Thus, descriptions of the 
FDTE and PFTG will not be repeated here. 

3. 3. 2.1 Error Analysis Processor - The EAP attempts to estimate 
the significant error sources of each of the shuttle nav sensors. 
The EAP follows the same design philosophy as the NFP . It reads 
the nav sensor data off an edited flight data tape. The EAP 
parameter base is made available through a parameter file tape; 
the contents of the file can be modified by the EAP under user 
control. EAP processing is directed through user- selected job 
control cards, controlling such options as processing start and 
stop time, IMU-to-be-processed selection, reference system 
selection, etc. 

3. 3. 2. 1.1 Source and Type of Inputs - Inputs to the EAP consist 
of the edited flight data tape (EFDT) , the parameter file tape 
(PFT) and a job control deck. The EFDT is generated by the flight 
data tape editor pre-processor. The PFT is generated by the 
parameter file tape generator pre-processor. The job control 
deck serves to specify user options and contains the following 
types of cards : 
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. function cards, specifying tape label and listable 
output page header information 

. data alteration cards, specifying modifications to 
the program data base contained in the parameter 
file 

. generation cards, specifying user- selected processing 
options 

. listable output control cards, specifying the type and 
frequency of program listable output 

. termination card, indicating end of the job control 
deck 

3. 3. 2. 1.2 Source and Type of Outputs - EAP outputs consist of 
tapes for further processing and listable output for visual 
analysis . 

All tapes are labeled with identifiers provided through 
the job control function cards. Tape output data consists of the 
state estimates and covariance matrix at each filter update. 

Listable output consists of a subset of the tape output 
data, provided at lower frequency. Such data might consist of 
position and velocity estimates and their computed variances, plus 
measurement residuals, in order to confirm filter convergence. 

In addition, error messages and diagnostics are provided as required. 

3. 3. 2.1. 3 Error Analysis Processing - A simplified flow diagram 
for the error analysis processor (EAP) is presented in Figure 3-7. 

The parameter file tape (PFT) is first read to set up the initial 
program data base. The processing of job control cards, which 
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Figure 3-7 Error Analysis Processor 
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generate header labels, set processing control flags and modify 
values in the data base, then proceeds similarly as for the NFP 
(see Section 3.2.2). If errors have been encountered in 
processing the job control card deck, the program terminates. 

If a new PFT is desired, the contents of the parameter 
file, as modified by the input data alteration control cards, are 
output to the specified tape. The data base parameters are then 
listed, following which they are converted from their colloquial 
type units used for I/O, to a consistent internal processing set 
( °/hr to rads/sec, min to rads, g's to ft/sec 2 , etc.). 

The EFDT is then read until the desired processing start 
time is found. The filter state and covariance matrix are then 
initialized. Also to be initialized is a set of INS data to be 
used as the basis for interpolating across any gaps which may be 
present in the recorded data. 

The program then enters its main loop. If more data is 
available in the EFDT, the next record is read. If the selected 
processing stop time has been reached, the program terminates. 

The filter is to compute the state transition matrix through 
the approximation 0=1+ FAt, where F is the system dynamics 
matrix, and At is the filter time step. In order to validly 
neglect the higher- order terms in the approximation, the filter 
time step must be constrained. Thus, if a data gap is encountered 
which exceeds the filter time step limit, a set of "pseudo-INS" 
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data must be generated by interpolation for use in the filter 
time transition. This interpolation would be of higher accuracy 
than the linear extrapolation used for the state transition 
matrix computation. The data to be used in this interpolation 
routine is kept current by updating it each time a record is 
read. 

Some elements of the state transition matrix are functions 
of integrals of specific force or integrals of products of 
specific forces (i.e., elements relating gyro g sensitive drifts 
to platform tilts). Given the rapidly varying character of the 
specific force vector for the shuttle, the allowable filter time 
step would be severely constrained. In order to avoid this, 
the specific force-dependent terms are accumulated at the rate 
at which the data becomes available, thus still constituting a 
faithful representation of the shuttle I MU specific force history 
when they are used at the more leisurely filter update rate. In 
this fashion, processing time is saved considerably. 

Next, if a measurement is available, or if the filter 
time step constraint is about to be violated, the filter state 
and covariance matrix are extrapolated to the new filter time. 

The order in which the filter is to process the measurement 
sensor data must be established. One such ordering, organized 
roughly from the least to the most accurate measurement type, is 
presented in the flow chart: barometric altimeter altitude, radar 
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altimeter altitude, TACAN VOR meas, TACAN range meas, MLS eleva- 
tion, MLS azimuth and MLS range. Such an ordering tends to 
process last the three measurements which are most sensitive to 
non-linear effects, thus using for them the best available state 
estimates. Actual processing of the measurements is, of course, 
dependent on whether the sensor data is available at the selected 
new filter time. 

For each sensor, a measurement correction and calibration 
step is first executed. This might include compensating the 
barometric altitude for non-standard temperature effects if 
temperature readings are available, correcting the radar altimeter, 
TACAN and MLS measurements for the fact that the sensor antenna 
is not co- located with the IMU, compensating elevation angle and 
range measurements for atmospheric refraction effects, etc. 

Following this, the appropriate measurement residual, measurement 
gradient vector and variance of the assumed measurement random 
error is computed. 

In order to avoid matrix inversions, the sensor data are 

treated as sequential scalar measurements. The following processing 

is common then to all measurements and is called for each measure- 

T 

ment. The filter estimated measurement variance (h P h + r, 
where h = measurement gradient vector, P = covariance matrix, and 
r = variance of measurement random error) is computed and used to 
apply a data-reasonableness test to the measurement residual. If 
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the square of the measurement residual exceeds some multiple 
of the estimated measurement variance, the measurement is neglected. 
Such a level might be 16 variances (or 4a) which contains over 
99% of the possible samples for a gaussian distributed random 
error. If the filter has been properly formulated, the covariance 
matrix should contain an actual representation of the la level of . 
the error in the estimates , and the test should therefore be able 
to edit out bad measurements. 

If the data validity test is passed, the measurement is 
used to update the filter state and covariance matrix. Depending 
on the measurement that has just been processed, and consistent 
with the pre-established processing order, the next measurement 
is selected, and the loop repeats until all measurements have been 
exhausted. Appropriate printouts and outputs to tape are made, . 
and the program loops back and reads a new EFDT record. 

When all EFDT records have been exhausted, or when the 
desired processing stop time is reached, the program terminates. 

3.3.3 EIP Summary 

The concept and formulation of the error isolation processor 
described in the previous sections provide the capability of 
determining magnitudes of particular navigation error sources for 
any one test. This processor is considered to be the least critical: 
it provides data that is valuable but not mandatory if the test 
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goal is only to determine overall navigation system performance. 
If additional detail is desired on subsystem performance, the 
EIP provides a method of obtaining this data. 
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4. POST PROCESSING DEVELOPMENT PLAN 


The development of the navigation data post processor should 
be divided into three phases : 

a. Specification Development Phase - where the preliminary 
design of the processor will be achieved; 

b. Im plementation Phase - where the design will be completed, 

and the software built and verified; and 

c. Test Phase - where integration testing will be achieved 
using data from the orbiter ground checkout facility. 

The rationale behind this plan structure and its relationship to 
other shuttle activity is presented below. 

pigare 4 — 1 shows the three phases and their relationships 
to the standard NASA development scenario and to the significant 
shuttle program milestones. The standard scenario is presented 
in the sections of the Figure 4-1 labelled "Development Activity 
and Documents" and "Review Milestones". The specification develop- 
ment phase requires an input of the navigation system requirements, 
design, and test requirements (as shown in the diagram) . With 
those inputs, it will be possible to achieve preliminary design 
review (PDR) of the post processor within six months. The design 
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requirements review (DRR) would be achieved between the second 
and third month, with the program plan developed in the first 
month. (The writing of the program plan could be taken as a 
separate project phase, if so desired.) 

The major reason for having a specification development 
phase separate from an implementation phase is the difficulty 
in being definitive in the scheduling and the costing of the 
software package implementation without a prior analysis of the 
navigation specifications. Hence, the program plan will be 
updated during the specification writing phase for the purposes 
of the implementation phase. 

The reason for a separate test phase (as opposed to including 
it in the implmentation phase) is the different phase logistics. 
Integration testing should be accomplished at the operational 
test facility using data from the actual vehicle ground checkout; 
whereas the other activities can be achieved at any location. 

A preliminary estimate of the level of effort for the develop- 
ment of various combinations of post-test processing tools is 
presented in Table 4-1. It should be noted that the level of 
effort required for the development of more than one form of post 
processor is not the sum of the man-months required for the 
independent development of individual processors. This occurs 
due to the similarity of certain modules in the various post 
processors . 
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Table 4-1 


Level of Effort (Man-Months) for Combinations of Processor Development 


PROCESSOR (S) 

SPECIFICATION 

DEVELOPMENT 

IMPLEMENTATION 

PHASE 

TEST 

PHASE 

TOTAL 

QLP 

6 

32 

4 

42 

PFNP 

6 

32 

4 

42 

EIP 

6 

32 

4 

42 

QLP & PFNP 

ii 

56 

6 

. 73 

QLP, PFNP) 





& EIP j 

16 

- 80 

6 

102 


There are two shuttle program forcing elements that act 
on the post processing development plan (see the bottom section 
of Figure 4-1) . The primary forcing element is the orbiter 
vehicle ground checkout. At a point in time when navigation 
equipment has been installed, the post processor should be ready 
to demonstrate its capability. According to the schedule in 
Figure 4-1, 14 months prior to that installation should be the 
post processing project initiation date. But project initiation 
requires input of a set of navigation specifications — — the second 
forcing element. Orbiter ground checkout begins in March, 1976, 
with the first captive flight scheduled for February, 1977. The 
navigation system should, therefore, be installed by the summer 
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(say, July) of 1976. This means that May, 1975, would be an 
appropriate starting date for the post processing development 
activity, with PDR occurring around November, 1975. 
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5. SUMMARY AND CONCLUSIONS 


The development of navigation data post-flight processors 
can be, if properly planned and implemented, a significant cost- 
saving device for the overall shuttle program. Using flight 
data from the initial test flights, it is possible to accomplish 
significant design tradeoff analyses of other flights through 
inexpensive simulation. Those simulations will have the added 
advantage of using actual operational data. The cost savings of 
simulation over actual flights are obtained, while the major 
argument against simulation is answered through the use of actual 
flight data. 

To summarize, the post-flight processing project outlined 
in this report achieves the following capabilities: 

. Rapid reporting of flight performance 

. Central control, synchronization, and standardization 
of data 

. Evaluation of alternative navigation schemes and modes 

. Performance accuracy analysis 

with the following cost and performance advantages to the overall 
program; 
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. Early identification of flight anomalies, hence a 
rapid initiation of remedial design allowing for 
tighter schedules between flights 

. Optimization of flight navigation configuration 

. minimum software requirements 

. minimum calibration and alignment requirements 
. minimum hardware requirements 

. Checkout of software without complex flight tests 

. Timely response to data requests from analysts, 
designers, and managers 

. Tight control and catalog of data inventory 

, minimum post- flight processing activity 
. no loss or damage of data 

. Standard programming language for all analyses 

. Common data base for all users 


The anticipated development schedule for this project 
(see Section 4) should be begun by May, 1975, in order to meet 
the overall schedule objectives of the shuttle program. 
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These equations can be readily derived from the fundamental /I 3 
equation of dynamics for a free material point written in vector 
form 


mW = F, 


where W is the acceleration of the point; and F is the resultant 
defined by the following sum: 


Here f is the resultant force of attraction of the celestial body; 
itlFi is the auxiliary mass force caused by the influence of the 
noncentrality of the planetary force field. This force must be 
taken into account for motion near the planet; and F 2 are the 
forces that do not have a potential function, for example, the 
force of air drag, which is taken into account for motion in the 
atmosphere . 

After canceling out m, the equations become 

A’] m 

j Z\ ' 

= ■ F 1 V rtij ALTT - iL + F w + ■ (1*1) 

H m i 

jmA A' m 


The system of differential equations( ( 1 . 1) of sixth order is 
thermo st general system of equations of motions of artificial 
spacecraft, including artificial Earth satellites. For assigned 


3 



initial conditions of motion (x 0 , yQ s zq ) , it must be integrated 
by numerical methods. In a special, but very important case when 
the point travels in the gravity field_of only one planet (j) 
and the perturbing forces F^ and (l/m)F2 are not taken into 
account, the system of equations (1.1) can be integrated in the 
final form. 

Let us transfer the origin of coordinates to the center of 
attraction, that is, let us jet xj = yj = Zj ~ 0 (we can neglect 
the influence of the mass of the AES on planetary motion), then /l4 

equations (1.1) take on the form of equations of unperturbed 
motion: 


x-\-k 2 m f -^- = b; 

3 r 3 

y + -~=0; 


where r - /x2 + y2 + z2" is the radius-vector of the point under 
study relative to the center of attraction. 

Let us multiply the first equation by y, and the second by x, 
and let us subtract one from the other. Then, 


■*d dy 
xy — zx~ — I x -Z- 

di \ dt 



By integrating this equation, we get 


In similar fashion, 
the first and third 


dy 

dr 


dz 

dt 

dx 

dt 


1 

il 

P 

(1.3) 

second and third equations, 
, we get 

and from 

, d ‘y —r • 

--3 ^ -C„ 

(1.4) 

1 

1! 

P 

(1.5) 


These three first integrals of system (1,2). are the integrals of 
areas . 
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If we multiply integrals (1.3) -(1.5) by z, x, and y, respec- 
tively, and if we add the results, we get 


C x x -\-Cyjy -\-C z z — 0, 


(1.5’) 


that Is, the motion of a material point acted on by .the central 
force applied at 0 occurs in a plane passing through point 0. 

This is physically clear since the side forces with respect to 
the, plane containing the . radius-vector r ; 'of the material point 
and Its velocity vector r are absolute. The position of this 
plane in space is entirely defined by the initial conditions of 
motion, that is, by the initial coordinates of a point and by 
the initial velocity of the coordinates with integrals (1.3)-(1.5). 

By moving in the plane (1.5?) > the material point preserves 
its;;sectorial velocity constant. We can easily see this since /15 
the left-hand sides of Eqs . (1.3)-(1.5) are the projections of 

the vector r x r onto the coordinate axes. According to the 
definition, the vector product of two vectors is equal In magni- 
tude (or modulus) to the product of their moduli by the sine of 
the angld between them, that is, equal to the area of the 
parallelogram constructed on these vectors. 

. In our case, r is the radius-vector of the material point 

arid r is its velocity. The area of a parallelogram constructed 
on these vectors is numerically equal to double the sector 
velocity. This means that under the effect of the gravitational 
force of one center (one planet), the material point will move 
along a plane curve, preserving its sectorial velocity constant 
and, therefore, preserving its projections onto the coordinate 
axes constant . 

Let us denote the sectorial velocity, that is, the increment 
in area A swept out by the radius-vector of a traveling point 
per unit of time by dA/dt. Let us find the value of the constant 
1/2 C, which it is equal to: 


where 




dA __ 1 r 
dt 2 

~ ^y 0 ? +(Vo - Vo ? + ( x oy<i - y.0 

c= j/q+q+q. 


( 1 . 6 ) 
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After integration* we get 


that is* the area of sector A increases in proportion to time t 
(to is the initial moment of time reference). Obviously, the 
projection of the area of sector A onto any of -“the coordinate 
planes will vary according to the same linear law. The constants 
C]_, C 2 » and Co are the projections of the doubled sectorial 
velocity C onto the coordinates planes Oyz, Oxz, and Ozy. If C^, 
C 2 , and C 3 are known, then therefore not only the sectorial) 
velocity 1/2 C is known, but also defined is theoorientation of 
the plane in which the point (of the AES) will move. The con- 
stants Ci, -C 2 , and C 3 are usually replaced by the more graphic 
parameters C, i, and ft, of which: i is the inclination of the 

orbital plane to the principal coordinate plane Oxy (if we con- 
sider the equatorial coordinate system, then the inclinationoof 
the orbital plane to the equatorial plane); and ft is the angle in 
the plane Oxy between the axis Os and the line of intersection' >b'f 
the planes of the equator and the orbit (nodal line), where the 
point at which the satellite travels from bottom upwards (as- 
cending node) is selected for the reference of the angle on the 
nodal line. 



Pig. 1.1. Scheme 
of motion of AES 
relative to se-, r 
lected reference 
system. 


From the spherical triangles (see /I 6 

Pig. 1.1) xnK and ynK (the arc nK is 
equal to 90 °) we find 

cos a=.-sin ftsin /; < 
cos —cos 9 sin i. 


Therefore , the integrals of areas (1 .‘3 ) - 
( 1 . 5 ) with new constants will be of the 
form 


yz — zy~C sin ft sin i; 
zx — xz ~ ~C cosft sin /; 
xy—yx—C cos/. 


Key : a . Orbit ; 

b. Equator; c. AES The longitude ft of the ascending node will 

vary within the limits from 0 to 3 6 0 ° f, and 
the orbital inclination i will vary from 
0 to l80°w If 0 < i < 90° satellite motion 
will be rectilinear, that is, it will be executed in the same 
direction in which the longitude Increases; if 90° < i < 180° , 
the motion will be retrograde. 
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a 


Fig. 1.2. 
For deri- 
vation of 
law of 
areas . 


The system of equations (1.2) gives us yet 
another integral — the integral of energy. In 
order to obtain this integral, let us multiply the 
first equation by x, and the second by y, and the 
third by z, and let us add the results, resulting 
in the equation 

. *■*+£/# + zz-\- {xx+yy+zz)=0. ' 


Integrating, we arrive at the following equality: 


4 


k-m , 


\ x* + if- + * 2 




where H is the constant. 

The expression in the parentheses is v2 — the velocity 
squared; therefore. In the final form the integral of energy 
will be 


1/2 k^mj ^ 

2 r 


Along the orbit the sum of the kinetic and potential energies 
of the AES as It moves in the central field remain constant. By 
using the integrals of area and energy, we can solve the problem 
of the motion of the material point acted on by the central force 

[58]. 

Let us write out the law of areas (1.6) in the coordinates 
r and <j> (angle of rotation of the radius-vector r in the orbital 
plane ) : 


r 2 ?=C. 


(1.60 


The same result for these orbits is obtained directly from Fig. /17 
1.2, namely " — ^ 


r{r-\-± r )^L-. 

At 


C-=2A 


As At + 0, this expression is transformed into (1.6?), 
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0 • ? p p p 

Assuming - r + r <J> and k mj = y we can write the equa- 
tion of energy as: 

T+ 2 ¥ ; 

Canceling out <}> and dt by means of Eq. (1.6 T ) we get 



Let us introduce the variable z = C/r, then 


or 



Furthermore , we introduce the constant e : 


2// + = 

C2 C2 


The constant e enables us to reduce the equation obtained 
above to the form 



dz 

— 4- d'$. 

1 

/£--( 

-*)' , 

After.' integration 

we get 


?+£>“ f- 


z C 

— arc cos — - 


u \ 2 ^ _ 

J1 

c) c 


z — 


cos (tp-f-D) 


and 


Z — yr~yr e COS (<P + D), 


where D is the fifth arbitrary constant. 
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Setting p = C^/y 


and returning to the variable r, we get 


/18 


and 


J L — J-e gos ftp -\-D) 

r p p 

r — - . 

1 + e cos (<p -f O) 


(1.7) 


This is the equation of a second-order curve in polar 
coordinates r and , related to the focus. In it e is the 
■eccentricity and p is the parameter of the curve. 

For the ellipse we have 

p=a e (l-e 2 ),e<l,H<0; 

for the hyperbola 

p = a^Ce 2 - 1), e > 1, H > 0, 


where a e is the semi-major axis of the ellipse; and a^ is the 
actual semi-axis of the hyperbola. 

For the parabola e = 1 and H = 0. 

When the angle <p = = -D, the polar radius r will be at a 

minimum. This is the perigee of the orbit if the Earth is the 
center of attraction . When <j> = tj) a - -D + tt; the polar axis in 
the case of an ellipse will have a maximum. This is the apogee 
of the orbit. 


The. angle = <j> + D = 4> - cf)^ is called the true anomaly. 

It remains to find the dependence of 0| on time. It will 
differ for the ellipse, the parabola, and the hyperbola. 

Based on the equation 


where 


r z dy=Cdt, 
_ p 

1 4- e cos ft 


C — 


(1-7 1 ) 


we get 
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db 

(1 4- e cos 8)2 


V 


( 1 . 8 ) 


„3/2 


dt. 


Instead of we Introduce the new variable a = tan JL-j 1 ; 
then -2| 


0— 2arctgri; . 




2 df\ 

tt? 


cos ft= 


l--n 2 

1 4* T) 2 ‘ 


The.; left-hand side of Eq. (1.8) can be transformed thusly 


where 


db 

( 1 4* e cos 8)2 


(1 + 


vo( 


2df} 

1 4- e 


n2\2 


1 -j- f)5 


2(1 + i\2)d^ 


2(14- t] 2 ) dt\ 

[1 + y] 2 + e( l_V)J 2 

1 4“ 


[1 + e + (1 — e)tfp (i 4- e) 2 (1 + Y7]2 )2 

Y=— - 
1 4- e 


dr\ t 
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With reference to this transformation, the initial equation 
(1.8) takes on the form 


1+7,2 dr\ ^< 1 + g > 2 
(1 4- yV) 2 2 fM* 



(1.9) 


Depending on the magnitude and sign of y, ,the solutions (or 
integrals) of this equation will differ. The most practical value 
is found in the case y ^ 0, which breaks down into two variants: 


and 


Y<0, e>\ ! 

y>0, e<\ . i 


However, it is simplest to integrate Eq . (1.9) when y - 0 (parabo- 
la) . For it , we get 
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Integration yields a formula for calculating the time of flight 
along a parabola 



{t — t) = tg ~ \~ 



i 

i 

l 

i 


where t is the time of transit of the spacecraft across the 
orbital perigee ( ft] = 0) 


For the ellipse and the hyperbola, the derivation of the 
dependence of the angle $| on time t by means of variable n is 
both complex and is an artificial process. However, referring 
to the fundamental nature of the resulting solutions, we present 
it in full [ 58 ]. First of all we note that 


there fore 


1 + *i 2 
(i -K vn 2 ) 2 


•fC + fl’i + i-j 

(1 + Yl 2 ) 2 


1 1 | V-l 1 1 

Y 1 -4- vn 2 Y (1 + YV) 2 ’; 


C (i + __ 1 r 

J ( 1 + Y 7 ! 2 ) 2 V J 


1 Y— 1 C 

1 *+* Y 7 ] 2 Y J (1 + Y 7 ! 2 ) 2 * 


( 1 . 10 ) 


Let us apply to the first integral in the right-hand side of /20 
the equation the formula of integration by parts 


hence 


dn __ T 
1 + y t i 2 1 + vn 2 

s + 2 \-i 

tr X 


A 




+ W 


(1 + Y 7 } 2 ) 2 
& 1 


] -I- y 7 ) 2 


2 1 


(t 4- Y 7 ! 2 ) 2 


(it j 

(1 + Y 7 ! 2 ) 2 " 


1 + Y 7 ) 2 




*1 

I 4- Y 7 ] 2 


Then Eq. (1.10) becomes 

[ ( ] + T l^) rf7 l __ V — 1 Tj 1 Y 4- I r dj 

J ( 1 4- Y 7 ] 2 ) 2 2y 1 4- Y r ( 2 2y J 1 4- y 7 ) 2 


( 1 . 11 ) 
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Using (1.9) and (1.11), we can write the integral of Eq. (1.8) 
as : 




1 Y-l 


n3/2 


(1 + e) 2 y 14 Y 7 ! 


Y 4 1 j 1 rfg 
2 ' y • 1 + Y 7 ) 2 


Since 


’ Y 


1 — 2e , 

“U ^TTT ' *> 


Y + 1 


1 — e 


then (1.12) can be written thusly: 


v7 

*3/2 


(/ t ) : 




/r»“ (1 + «)(!- * 3 ) I J l+YV 1 + Y 7 ] 2 


rf-»l 


— <? 


( 1 . 12 ) 


(1.13) 


From Eq. (1.13) we directly get the formulas for calculating 
the time of flight in the case of an elliptical and a hyperbolic 
orbit. 

For the ellipse (e < 1, y > 0) , we have 


\ - = -4^arctg (ii /y). . 

J 1 4 Y 7 ! 2 / Y 

We introduce the new variable E, using the equality 


i/v-tB-f-. ; u.i4) 


Therefore, we have 


£*=2 arctg (ti )Ay); 


1 -f- Y 7 ) 2 


2 

/y 


, £ £ 1 

sin — cos — = sin 

2 2 / Y 


E 


and Eq. (1.13) becomes 


v v- U — T ) 
a f( l~^ 2 


2 

(1 I <.')(!* -* 2 ) 


|— r^arctgr} /y ^r-sinf 

IKY 2 > v 


/21 
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or 


V H- 
nW 


,{t — x)—E~ *. sin E, 


(i.13 1 ) 


where E is the eccentric anomaly. 

Eq. (1.13’) for calculating the time of flight in an 
elliptical orbit is called Kepler's equation. 

The physical meaning of E is clear from Fig. 1.3, in which, 
besides the trajectory of the satellite motion in an ellipse, is 
shown the trajectory of motion along a circle of point A that 
has an x-value that is identical with the satellite: 


FB=zOB — OF\ or 


r cos 6 1 — o e cos E — a^e. 


But from the equation of an ellipse it followstthat 


1 4- e cos 9 


Dividing the two last equations by each other, we get 


COS + gcos 8)(cos E—e) 

1 — <22 


Therefore , 


or 


cos E 


e -f- cos 9 
1 -f e cos 8 


e -r 


tg2 T 




+ tg2 T 


1 4- tg2 


B 


1 4- e- 


i-«g 2 v; 


1 4- tg‘ 2 ~ 


14-e — (1 — e)tg2 — 


1 4-<?4-(l~<?)ig 2 — 
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formula re- 
lating the 
angles of the 
eccentric ano- 
maly E and the 
true anomaly 

Key: 1. AES; 

2. Orbit; 

3 . 3-0 


After obvious transformations we have 


l + tg a -f-+e+etg*|--tgfU+ e tg*-l = 
= l-tg a |-+e-*tg a -|+^.±- e tg a -i; 
(l + e)tg a Y=(l-'e)tg a Jt 


/22 


and, finally. 



\ 

(I.l4«) 


Thus, we get Eq. (1.14 1 ), which gives us a 
relationship between E and . 

For the hyperbola (e > 1, y < 0), we have > 


1 y ^]2 


l In 1 + ^ ~~ yn 

2 / — Y 1 — ‘ 


We introduce the variable q, defined by the equality 


tg~ => — Y'H — 


y \ + e & 2 


After obvious transformations of (1.13), we get 

(/ - t)= e tg q - lntg (-5- + -§-)■ 


or 




t — x=-~{e$hq— q). 
Vv- 


Thus, the problem of integrating the equations of unperturbed 
motion (1.2) has been solved. As to be expected, its;:solution 
depends on . six parameters: a, e, T,:;'.i, fl, ^andj.^ (or w), and in 

rectangular coordinates — on xq, yg, zq, xq, yo , and zn. The 
trajectories are in the form of ellipses (periodic solution) or 
hyperbolas (aperiodic solution), A singular case is the parabola 
(e = 1) and a particular case of the ellipse is the circle (e = 0) . 
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For the circular orbit we have 


:&=<p; a= R; *— 

/i* 


where <j) is the center of 'the angle. 

The period of revolution for elliptical orbits, as follows /23 
from Kepler's equation is 


r=-^=-a 3/2 i 

Y v- 

and does not depend on eccentricity. 

We note that for elliptical orbits Eq. (1,8) can be readily 
integrated by means of the substitution 


, & / \ + e , 

tg T = y “ tg 2 • 

By differentiating this formula, we get 


sec 2 — = 

2 


or 


/ 1 +e 
' 1 —e 

2 £ 

OCL ___ 

2 

f \ + e | 
1 — e ' 

'i+tg 2 


Hence it follows that 


db=* 


Y\ — e2 

I — e cos E 


dE. 


On the other hand, from the equation *cos E = — & ; we have 

1 + e cos 0 ; 

1 l — e cos E 1 

1 + e cos $ l — d 1 

Using this. .expression, we can define the radius-vector of 
an elliptical orbit in terms of the eccentric anomaly, that is, 

r~a( 1 — -ecos E). 
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Further, from Eq. (1.8), with reference to these formulas, we 
have 


v't* 

aW 


dt—{ \ — e cos E) dE. 


Integration gives us Kepler's equation ( 1 . 13 T ) . 


Calculation Formulas for[ Unperturbed Elliptical Orbit 

Above we obtained the following fundamental formulas, which 
can be used ".in calculating elliptical motion: 


V- 

~2 


r 2 d = C * 

1 


integral of areas 
integral of energy 


ge’neral formula for any motion /2 4 
expressing the radius vector of 
a satellite in terms of p, e 
and ; 


3/2 _ is Kepler’s equation for cal- 

t~x — — — (E — esinf)' culating the time of flight in 
y ** elliptical motion; and 



is the parameter of the orbit 


a l/i , 2 HC* 

and e = y 1 i — 

city. ** 


is eccentric 


In practice it often becomes necessary to use a series of 
formulas deriving from the fundamental relations or obtained from 
geometrical considerations. We present some of them. 

We will assume that the point where the satellite is inserted 
into orbit has the coordinates: geocentric latitude and 

longitude Xq. The azimuth of the absolute satellite velocity at 
this point will be denoted by 6q. From these quantities we can 
calculate the inclination i of the orbit to the equatorial plane 
and the argument of the latitude at the moment of insertion, 
determining the distance of point B of the insertion into orbit 
from the orbital node (Fig. 1.4). This argument ( (uq)- is> expressed 
by the formula uo = to + (ftJo > where to is the argument of the perigee 
( 7r ) of the orbit. 1 
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Examining the spherical 
triangle AB C, we can obtain 
the following formulas: 


cos /= sin 8 0 cos ']> 0 ; 


sin h 0 = 


sin j 
sin i 1 


Ctg« 0 = 


cos 5(> 
tg^O ’ 


(1.15) 

( 1 . 16 ) 


Pig. 1.4. Elements of orbit of 
an AES . 

Key : 1 . Orbit ; 2 . Direction 

toward perigee (tt); 3 . Orbital 
node ; 4 . Equator 


The angular distance of 
meridian of the point of 
insertion from the orbital 
node can be defined from the 
same triangle: 


sin vq — sin o 0 sin « 0 ; 


cos v 0 = 


COS Uf) 
COS 'J'o 


(1.17) 

( 1 . 18 ) 


We will measure the longitude of the orbital node from the 
point of vernal equinox. Here it can be determined as the direct 
ascension of the orbital node. At the moment of insertion into 
orbit, we have 


' -0 + Vo — S 0 * 


(1.19) 


where sq is the sidereal time at the meridian of the point of 
insertion B into orbit at the moment of insertion; and vq is the 
arc from the ascending node of the orbit to the meridian of the 
point B. 

Hence we can determine the arc Aq. Obviously, in order to 
obtain the longitude of the node X measured from some terrestrial 
meridian, it Is necessary to add to ftp the quantity (~sq t ), where 
Sq t is the corresponding sidereal time for the selected meridian. 

The resulting formulas (1.15-1.19) are valid, of course, not 
only for the point of Insertion, but for any point on any orbit. 

Now let us find expressions for the absolute coordinates 
x, y, and z, and their time derivatives in terms of the orbital 
elements ft, to, i, a, e, and t* To do this, let us consider 


the spherical triangles xAB, yAB and zAB formed by the coordinate 
axes and the arc AB. Using the theorem of sines for the direction 
cosines of the radius-vector r of point B, we get 

cos (rx ) — cos it cos 2 — sin u sin 2 cos i; ' 

cos cos u sin 2 sin it cos 2 cos /; (1.20) 

cos {rz)- ~ sin u sin i. 


Therefore , 


AT--r(cos u cos 2 ~ sin u sin 2 cos/); 
y^=r{cosu sin 2 -(-sin «cos 2 cos/); 
z-~r sin « sin i. 


( 1 . 21 ) 


Differentiating these equations with respect to t, we have 

x=r~ — r(cos 2 sin w-4-sin 2 cos ttcos/) 0 ; j 
r ' j 

i y—r— — r ( sin 2 sin u — cos 2 cos u cos /) 9; \ 

\ r [ 

I z~r cos u sin /4>4- — r. 

( 1 . 22 ) 


Here 




The resulting formulas (1.21), (1.22), (1.13*), (1.7') a and /26 
and (1.14), that is, the expressions of the absolute coordinates 
x, y, z, and their derivatives x, y, and z in terms of the:'- 
orbital parameters (elements) 2, w, i, a, e, ,x..dnd time t are the 
solution to the system of equations of unperturbed motion of a 
cosmic body (Earth satellite). For each point in space at which 
the satellite is located, from formulas (1.21) and (1.22), by 
knowing the five geometric parameters defining the orbit and the 
one kinematic parameter (x),for any moment of time we can calculate 
the absolute coordinates and the satellite velocity. 

Note that often instead of parameter x the more convenient 
parameter :.M o' appears. 


Let us denote 


M = E - e s i n E = n (t _ t v ) + M „ , 



where 


M is the mean anomaly; 
n^Y^a-w is the mean motion; 

to is any specific (initial) 
instant of time; 

A/ 0 =/z (/"o — T ) is the mean anomaly for the 

instant tg 3 the so-called mean 
anomaly of the epoch; and 

Mq is the parameter replacing 

t and characterizing the value 
of M for the instant tg- 

Let us find the expressions for the orbital parameters in 
terms of the elements at the end of the powered flight section: 
vo a 9o (angle of inclination of velocity vector to the horizon), 
and r 0 . 

Based on the formulas 


r= r 2» — C; p = — 

1 -i- e cos & p- 

the complements of the velocity vector can be expressed thusly: 
in terms of the radius-vector 



— e sin 9 
P 


in terms of the transversal 

<v n = — — (1-j-^cos ft). 

Hence the modulus of velocity is 


( 1 . 23 ) 


( 1 . 24 ) 


^ 2 =“^n + 2£cos&+e a ) ; 


( 1 . 25 ) 
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/27 


the velocity at the perigee (maximum) is 

(1.26) 

and the velocity at the apogee (minimum) is 




(1.27) 


For the radius-vector of a satellite at the perigee (mini- 
mum) , we have 


r K —a{\—e)\ 


( 128 ) 


at the apogee, we have 


i 

r«— «(! + «)• ‘ 


(1.29) 


Substituting the 
with reference to the 


values of v^ and r^ into the energy equation 
fact that p = a(l - e 2 ) , we get 


Jh. (1.30) 

2 r 2a , 

Thus, total satellite energy depends only on the semi-axis a of 
the satellite orbit. 

Replacing in Eq. (130) v and r by vq and r 0 at the point of 
insertion into orbit, we get 


where 



(1.31) 


The quantity k is the square of the ratio of velocity vq at 
the initial point to the velocity in the circular orbit with 
radius rp . 
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' 2 

Actually , when 1 , , based on Eq . (1.31) a = r 0 , 

which corresponds to the circular orbit. In this case. 


*~i/ t = l/T =const - 


The flight velocity v c j_ [ci = circular] in the circular orbit 
(r = const) is constant and can be expressed by the formula: 


Vci =j/f. 

This velocity is called the circular velocity at the dis- 
tance r from the center of attraction. 

If in Eq. (1.31) we set k = 2, the semi-axis of the ellipse /28 
a becomes equal to infinity, that is, the elliptical orbit is 
transformed into a parabolic orbit. Here we have 


v par 


which is called the parabolic velocity at distance r from the 
center of attraction. 

If we take r = R (where R is the radius of the Earth) in 
the formula for the circular velocity, we get the first escape 
velocity 


and on analogy 


is called the second escape velocity. 

Let us find the expressions for the orbital element p and e 
in terms of the initial parameters ro 3 v 0 and 0 q defining the 
ellipse. Since 
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r^Vr.— V pp = rv cos 0 


(where 6 Is the angle of inclination of the velocity vector to 
the local horizon), the orbital parameter is 


p=r 0 k cos%. (1.32) 

Replacing in Eq. (1.32) parameter p by its expression 
r— ~(1 , we get the following formula for eccentricity e: 

l R | 


e=Y(k— If cos a 0 o -fsin 2 G o . 


( 1 * 33 ) 


The angular distance (true anomaly) of the point of insertion 
from the orbital perigee will be determined from the polar equa- 
tion of the ellipse by 


co$& 0 =— (&cos 2 0 o — 1), ' (1-3*0 

e 

and the angular distance of the perigee from the orbital node 
will be 


ih—Uq — 



( 1 * 35 ) 


If at the end of the powered section 0 O “ 0, the orbital perigee 
coincides with the end of the section, and here O 0 =0. | 

The apogee distance and the eccentricity are determined by /29 
the formula 




2 - k 

e---k- 1. 


The period of revolution of the satellite 

T L- ( r ° 'f 2 

i i * | p, \2 — k j 

depends only on the orbital semi-axis a. 


(1.36) 

( 1 . 37 ) 


( 1 . 38 ) 
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Selection of Position of Orbital Plane 


The position of the orbital plane relative to the Earth is 
defined by the latitude and longitude of the point of insertion 
into orbit and by the azimuth of the velocity vector at this 
point . 

The position of the orbital plane relative to the ecliptic 
and* therefore, relative to stars and the Sun is characterized 
by the inclination of the orbit to the equatorial plane and the 
longitude of node £2 q . The inclination of the orbit is uniquely 
defined by the azimuth ( 6 o ) of the velocity vector at the point 
of insertion and by the latitude (tpo) of this point. The longi- 
tude of the node is defined by these same quantities, and, in 
addition, by the time of insertion into orbit. The time of 
insertion must be taken into account, since the launch is made 
with Earth rotating in absolute space. 

For convenience, let us use sidereal time s which gives 
directly the angle of rotation of the globe relative to fixed 
coordinate axes. 

Thus, the position of the< orbit at any instant of time rela- 
tive to Earth and Sun is defined by the quantities Aq* <$0 

and sq. These quantities are selected based on the purpose of 
the satellite and its design features, and also several other 
factors. For example, if the satellite has solar batteries, then 
it is useful to specifically orient its orbital plane relative to 
the Sun. If the satellite has a celestial orientation system, 
these quantities should be selected with reference to the capa- 
bilities of its normal functioning. 

Let us dwell in more detail on selection of the orbit relative 
to Earth, which is significant for observation of the satellite, 
setting up communications with it, and for several other purposes. 

Depending on the parameters of the insertion into orbit, the 
rotation of the Earth has different effects on the initial ab- 
solute satellite velocity. Since at any point during the powered 
section the influence of the Earth’s surface velocity can be 
assumed constant, the increment in velocity in converting from the /30 
relative to the absolute coordinate system is 


Av m 0.465 cos sin 8 0 cos 0 O ^=0.465 cos i cos 0 O km/s e c , 


where the numerical coefficient is the Earth’s surface velocity 
at the equator. 
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Pig. 1.5. Selecting the posi- 
tion of the orbital plane. 


Key: 1. Direction toward peri- 

gee; 2. Equatorial plane; 

3. Orbital plane; 4. Orbital 
node 


The formula shows that 
in order to obtain the 
most distant orbit it is 
desirable that the latitude 
of the launch point be 
selected as close as possible 
to the equator, and the 
azimuth <5 q — as close as 
possible to 90° , 

Denoting by and Aj 
the latitude and longitude 
of point E on Earth, res- 
pectively, above which at a 
given instant tj the satellite 
will pass (the satellite is at 
the zenith over this point 
at instant t j ) , by inspection 
of the triangle DBE (Fig. 1.5) 
we can obtain the following 
formulas : 


sin A Jj=s sin % cos -f cos 4> 0 sin cos . 

(1.39) 

Xj — Aq 4“ A^y — A Ay ro t 

(1.40) 

cosaX. — c °s ?/ — sin +/ sin +0 . 

1 COS <\>j COS 4^0 1 

1 \ 

• 

1 — t 

sin Py sin B 0 

. sin *h= co»"^ — 

(1.42) 


where AA-? rot = wg(tj - to) and wg is the absolute angular velocity 
of rotation of the Earth [rot = rotation; E = Earth]. 

If the satellite, orbit is near-circular, it can be assumed, 
approximat e ly , that 



(1.43) 


When using Eqs. ( 1 . 39 ) - ( 1 . 43 ) , we must consider that if the 
azimuth of the velocity vector relative to a rotating Earth (6 q t ) /31 

is the initial azimuth, then we must convert to the absolute 
azimuth <$q by the formula 


f 


0,465 cos ^,i cus 
v cos 0o 


(1.44) 
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Using these formulas, the projections of the orbit onto the 
Earth were calculated for different azimuths at the point of 
insertion, whose first orbits, according to the change in angle 
3j from 0 to 360° , are shown in Fig. 1.6. In order to obtain the 
projections of the other orbits, we must shift the curves along 
the latitude westerly, by the quantity (for the lower orbits) 


15,22 [h-\)T\ 

T_ 

where n is the ordinal number of the orbit, and T n is the period 
in hours (h) . 


1.2. General Problem of Inserting Spacecraft into Orbit 

Let us examine the case of the insertion into orbit of a 
satellite by means of some two-stage launch vehicle (LV), whose 
first stage uses a ramjet engine (RJE) and whose second stage 
uses a conventional liquid-propellant rocket engine (LPRE) . 

The theory presented below can be easily extended to the case of 
a large number of launch vehicle stages. 

We know that the use of an RJE in the atmospheric section of 
a flight makes it possible to substantially improve the power 
capabilities of the LV through use of air-oxygen in the Earth's 
atmosphere as oxidizer. 

We will not dwell on methods of calculating RJE, advantages 
and disadvantages of RJE compared to other types of engines, 
and the difficulties confronting the builders of this kind of 
engine. We refer the reader interested in these matters to 
certain published studies (for example, to [45, 57]). 

Our problem will be determining the optimal law of control of 
the LV engine with RJE over the powered section, that is, deter- 
mining the optimal control program with respect to the angle of 
attack a and the mass flow rate of fuel B. 

Due to the complexity of the problem we are to solve, it is 
not possible to obtain analytic functions. Therefore, the method 
of optimization presented here has been developed as applied to 
electronic computers. 


Fdrmulation of the Problem /33 

We will consider the motion of an' LV under the following, 
assumptions and presuppositions. 

1. The launch vehicle is regarded as a material point. 
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Fig. 1.6. Projections of first orbits for different 
azimuths 60 at point of insertion. 
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uo 
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2. The motion of the LV is considered in the vertical plane 
coinciding with the plane of a great circle of the Earth. 

3. The amount of fuel in the first stage, and the total 
weight of the first and second stages, are given. 

4. The field of gravity is central, and the Earth is not 
rotating . 

5. The limiting values of the variation in the angle of 
attack a m ^ n and a max are small, as the result of which we will 
assume that cos a ~ 1, sin a = a (in the extra-atmospheric section 
the selection of the boundaries a m ±n and « max will be made by 
using the terms P cos a and P sin a instead of P and P a in the 
equations of motion (1.45)). 

With reference to the foregoing assumptions and presupposi- 
tions, the system of equations of motion will become 

V = i sin 8 = /,; 

m r? 

■ £_ cose + J^2LL =/l ; 

mV rW r 

£ = f2-l/cos9==/ 3 ; I (1.45) 

hz=V sin 0=/ 4 ; 

m—— / 6 , j 

where 

<?x==c* # (M) + cJ(M)a a ; M)a. 

Let the following constraints be assigned: 


g 1 = (a-a mtn Ka-a m J<0; k ; 

. (i-W 

C4 = ?-‘7max<0; (1.^9) 

e,=v-/(A)<0. (-1.50) 


Here and above a min , a ma x 3 Sminj Bmax, n ylm i n , n ylmaX3 and 
qjnax^ are some preassigned constant numbers, and £(h) is some 
function of altitude, assigned in advance on the condition that 
for motion along the curve V = &(h), the equilibrium temperature 
of the LV * s surface at some characteristic points^on the surface /34 
is equal to the allowable maximum value; 
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V is flight velocity 

0 is the angle of inclination of the velocity vector to 
the local horizon 
L is flight distance 
h is flight altitude 
m is instantaneous mass of the LV 
P is engine thrust 
a is angle of attack 
3 is mass flow rate of reaction mass 
S is characteristic area 
q is velocity head 0 
y - 3 • 98 • 10 5 km3/sec^ 
r is radius-vector 

rg - 6381 km — the mean radius of the Earth 
c x is the drag coefficient 
cy is the lift coefficient 
ny-^ is the transverse load 

M is Mach number 
P S p is specific thrust. 

The inequality (1.46) is a constraint on the angle of attack 
a, and the inequality (1.47) is a constraint oh the mass flow rate 
of fuel. Since we will take the angle of attack and the mass flow 
rate of fuel as control function, inequalities (1.46) and (1.47) 
can be called constraints on the control functions. Inequality 
(1.48) is a constraint on the transverse g-load n yi , which as we 
know is defined by the formula 




c*ct Sqr* 


As we can easily see, ny-j_ depends explicitly on the angle of 

attack a. However ny-^ depends also on the phase coordinates (for 

example r, m, and so on). Usually this kind of constraint is 
called constraints on the zero-order phase coordinates. 


Inequality (1.49) indicates that the velocity head q is 

M = 0 13. = 

3a u ’ 36 


constrained. We can show that 


0, while at the same 


time 


3a 


dq 


dt 


a constrain 


/ 0. Customarily, this kind of constraint is called 


on first-order phase coordinates . The inequality 
(1.50) is a constraint on the equilibrium temperature of the LV’s 
surface at some characteristic points thereon. It can shown that 


35 

3a 


5 _ 


9S 


= 0, 


3(3 


5=0,ai 

3 3a 


fdS, 


■dt 


f 6 .0. This indicates that constraint 


( 1 . 50 ) is also a constraint on first-order phase coordinates. 
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Let us take as the optimality criterion the minimum mass flow /35 
rate of fuel during the insertion section 


/ 



(1.51) 


Suppose it is required for the LV to pass from one stage at 
instant to 


V{t*)=V» 


0 ('o) ®oi ^o) — A>» A(/ 0 )=A e , m (^=s/7j 0 


(1.52) 


to the final state at instant tf.* 

V(T f ) = Vf, 0(t f ) = 0 fJ L(tf) = L f , h(tf ) = h f . (1.53) 


The problem is formulated as follows: among all the admis- 

sible controls a and 3 satisfying constraints (1.46) and (1.47) 
and converting the phase point from position (1.52) to position 
(1.53), we are to find those at which the trajectory V,,0, L, h, and 
m, corresponding to them and satisfying the differential relations 
(1.45) would give a minimum for functional (1.51), in order that at 
any point of the trajectory the inequalities (1.48), (1.49), and 
( 1 . 50 ) are satisfied. 

We will assume that the functions V, 0, L, and h are continu- 
ous and are piecewise-dif f erentiable ; the function m is continuous 
and piecewise-dif f erentiable , and that the moment of staging 
changes discontinuously by a specific quantity; the functions a 
and B belong to the class of piecewise-continuous functions; Bmln* 
Bmaxj a m in and a max are such that they are different for different 
stages . 

Since the LV we are considering is assumed to be of the two- 
stage itype, the thrust, aerodynamic, and other characteristics of 
each stage differ. This means that this problem is in the category 
of discontinuous variational: .problems of the second type. 

One feature of this kind of problem is that at the instant of 
staging the right-hand members of the system (1.45) suffer a 
discontinuity. We will assume that the surface of the discontinui- 
ty is specified by the equation 

f ~m — m 2 — \ (1.54) 


where m is the instantaneous mass of the LV; m 2 is the mass of the 
LV after staging of the first stage; and Amp is the mass of the LV 
expended in the separation of the first stage. 
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Method of Optimization 


To- solve the problem formulated above, let us apply the L.S. 
Pontryagin principle of the maximum [53]. Under the principle of 
theuimaximum, we introduce some function H: 


a -55) 

i-i 

where A]_ are certain undetermined multipliers. 

Then the first condition for the stationary character of this 736 
functional in vector form will become: 


z dxi 

(Xy—V, X Z 


oti +*-£-+*. -*H-* 


Ox i 


dxi 


oj 5 
Oxj 


0; x.,=L\ x, — h; x^m). 


(1.56) 


The multipliers ^3, ^4, and ^5 are defined from the equations 


dH 

du 


*1 


du 


"i " 



Oh 

du 


+' t ‘< 

(a— a, ?), 


££5 

Ou 


(1.57) 


9 W 

where £4 and £5 are the total derivatives in time of the functions 
£4 and £5. Eqs. (I.56) and (1.57) indicate the fact that the 
desired extremal in the general case is a piecewise-smooth curve 
consisting of a finite number of sections, each of which either 
lies at the boundary of a given region (1.46) — (1.50)> or else 
is within these regions. In the second case, the optimal control 
is defined from the condition 



Ou 


(1.58) 


In this case all the multipliers = 0 (j = 1, 2, 3, 4, 5). 

If the control obtained from Eqs. (1.58) does not satisfy 
certain constraints ( 1 . 46) -( 1 . 50) , it is determined from the condi- 
tion of belonging to these boundaries and the condition of a 
maximum of the function H. In this case, 9H/3u 7 0. This Indi- 
cates that one or two multipliers are nonzero (the maximum of 
multipliers simultaneously nonzero cannot exceed two, since 
there are two control functions), and a wholly determined multiplier 
ipj 7 0 occurs only when = 0 corresponding to it. For example, 
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the optimal control is a control corresponding ten. the boundary 
n Yl “ n y max , only the multiplier ^3 whose numerical value is 

determined from the solution of Eq. (1.57) is nonzero. 

According to Pontryagin's principle for- autonomous systems 
in free time tf, the exact upper bound to function H can take on 
the same value, equal to zero, along the entire trajectory, that 
is, at any instant of time the following equality must be satisfied: 

H" = H + (1.59) 


The sign. indicates that the function H is taken at the left, /37 

and "+" is taken at the right of the fixed point at the given 
instant of time. 


As for the multipliers X-^, at the moment the first-order phase 
coordinates arrive at the bounds of the admissible regions, and 
also at the instant of staging, they can suffer a discontinuity of 
the first kind. 


The conditions for a discontinuity of the multipliers X^ at 
the instant of insertion are as follows: 

— at the boundary ^ = 0 


V= x r+ v 


dxi 


( 1 . 60 ) 


— at the boundary ^5 - 0 


X+_ ^ 4- Vl 


0x ( 


(1.61) 


At the moment of staging, the conditions of discontinuity of 
the cofactors are of the form 


Xt — \ + V 2 


<?/ 

Oxi 


( 1 . 62 ) 


Conjugate System of Differential Equations 

As already indicated above, the characteristics of each stage 
differ. Because of this, the trajectory breaks down into individu- 
al stages separated by the surface (1.50), in each of which stages 
the function H is of the form 


H = 4* A5/2 4~ Xg / 3 4* X4/4 4“ A &./V 1 


(1.63) 
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Eqs. (I.56) in phase coordinates at each of these stages will be: 


i ■ '-l \l 
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^2 


5A, 

5V 

0/2 
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d/a 1 d/ 4 , 5$ 4 , iti 5g 5 

ir~ x *“ir +t4 ir + ^ir 


A„ 0, 
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Of i 
dh 
0 Jj_ 
dm 


dh 

d / 2 


dh 

5^3 


dU 


dh 
Ok 5 


5/t 


*+*'*+*' at- 
om Om dm dm 


( 1 * 64 ) - 


Here 


o/y 
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0V m l OV 
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50 

5/2 __ 
5A 

5/ 2 _ 
5m 

g/a_. 

5V ‘ 

5/3 . 

50 " 
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dh 
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OV 

5/4 
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dv 


=—(2n yt -ih, ml „- < V- V ' ) ; ' 

mu. 
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(2tt 4 - t ^ifrnilu n 


aa 
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mu 
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dV 

2 

dA 
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1 

dQ 

d0 

2 

dA 

d$4 __ 

1 

dSQ 


g/l 

dK 


do 


Oh 2 aA3 

f ^4 ^/l . 
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de 
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dK dP dA 

dSs .... d/] d/ (A) 

dfl dO dA 
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dA dA 

d$s _ d/i 
dm dm 

dc. r 1 ( dc 
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dc 


dM 
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As indicated earlier, the LV we are considering is-.! a two-stage 
LV^using a RJE in the first stage, and a LPRE in the second. 

As we know, the specific thrust of an RJE for a selected 
engine and specified fuel components depends basically on the actu- 
al excess air coefficient a<h and the Mach number M, that is, it 
can be approximated by a polynomial of the n-th degree 


^ I =^ l) (M)+^i(M)a (t) +^ 2 { M )a2 4-. . . + A„(M)aJ, (1 . 65 ) 
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where 


a o=- £ r 1 -- o 0 =^-;* 1 =a(M)a+*(M); 

P X 

Aq (M) , A]_ (M) , A 2 (M) , . .., A n (M), a(M) , b(M) are some functions of 
the M number; k]_ is the coefficient of flow compression in the 
LV body; F is the area of the diffuser inlet; 0 is the coefficient 
of diffuser mass flow; p is the atmospheric density; and k is 
the stoichiometric ratio . 

Let us limit ourselves to examining only the first two terms 
of series (1.65). Then the thrust P will be of the form 


P=Z»\A x D,{aa-\-b)+A 0 % 


( 1 . 66 ) 


The thrust of a LPRE, as we know, is determined by the formula /40 


P^=g 0 {J u '?~(.\>S a Q), 


(1.67 ) 


p 

where go = 9-81 m/sec is the equation of terrestrial gravity at 
sea level; Jo is the specific thrust of an LPRE (in vacuo); S a is 
the area of the nozzle discharge; and eo = RT/gQ (R - 287 . 05 m 2 / 
/sec 2 * deg is the gas constant per kg-mass of air, and T is air 
temperature ) . 

Then the partial derivatives in P will be: 


— over the section up to the separation of the first stage 


OP 

dV 




M( 


6H- 

1 0- f , 1 

\ 

a*f dy\ a*A\ 

t?M / 

/ da 

f a 

+ £)]H 

t. 

1 - ■■ U 

<JM r 


L-j(aa-(-6) + 




— over the section after separation of the first stage: 



1 

1 


dP 
Oh ~ 


dh 


In solving Eqs. (1.64), we everywhere assumed that the speed 
of sound a* is constant in altitude. It was also assumed that the 
functions p, l{ h), <j>, c a , c x q, a, b, Aq 3 and Ap, and also their 
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<?Q d 2 Q al (A) omji) dc 

partial derivatives oh ' Oh*' oh ’ ’ <>M ’ 77 - 


A -2., £di and «dL 

dM t>M f>M <?M 1 r)M 


are known in advance and specified in the form of tables or ap- 
propriate graphs ♦ 


Determination of Optimal Control 

At any fixed moment of time x, the phase coordinates x-j_ and 
multipliers (cofactors) A^ are fixed. Therefore the function 
(I. 63 ) at this instant depends only on the controls a and 3. We 
also need to find a and 3 such that the function H takes on its 
exact upper bound. 

We will show that the constraint on the transverse g-load n 
can be reduced to a constraint on the angle of attack a. ^1 

Actually, in order for the transverse g-load n y ^ at any in- /4l 

stant of time x not to exceed its specified value, it is necessary 
that the angle of attack a lie in the interval 


where 


® ^ ®1niax? 


a linln — 




cfSqr 2 ’ 


Imnx 


n. 

h xmax 


mpi 


c;s q r* 


Therefore , 
constraints are 


it can be assumed that at any fixed moment two 
imposed on the control function a: 


iln 


< a < a, 


• 1 mi 


n <a<a 1( 


The last two inequalities can be replaced by the single in- 
equality of the form 


where 



a 2min 

^ a (Jo max » 


a 2mln — 

= a min> 

if 

®min 

> 

a lmin> 

^2min ~~ 


if 

^mlii 

< 

^linin’ 

a 2max ~ 

“ a rmix* 

if 

a max 

< 

Q linaX’ 

^Sinax ” 

“ a iniax) 

if 

a max 

> 

a imax* 


( 1 . 68 ) 


(1.69) 
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The function H can be represented in the form: 


— over the section up to the separation of the first stage 

/ 

H=A u a? + E u a+E n ap+N u il + T n ; ' (1 - 70) 

— over the section after the separation of the first stage 

H = H i2 a2 + -f N^i -f T n . (1.7D 


Here 


• 4 ‘, = ir (-f ~ Kc * Sg } 

E n = ~ + cjS?); 

N n =* — g«A„-\\ 

m \ 

T =x sin o\ + — cos s( V M+ 

\ m r- j r \ rV J , 

I 

+ X 3 -^-Keo3 0-j-X 4 Vsin0; I 

I 

m " ' 

I 

Eiss c’Sg - g„<?oS 0 o ) ; 

Fa= iv evI “ ; 

^12 ”'T~ So^o ''■S' 
m 

T a = ->. i (&&M±l^ + JL S i„t)±±L CO so(v -■£-) + 

V m ■ ' r l rK j 1 

-f a 4 V sin S-j-XgH-l/ cos & . 


Similarly, let us represent the relations imposed on the 
trajectory for motion along the boundaries £4 = 0 and £5=0. For 
the LV to move along these boundaries, the following must obtain: 

Ok -and £5 = 0. 

In expanded form, this can be written as: 
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— over the section prior to separation of the first stage 


£ 4 — -f £’ 21 ° ^2ir j “ L > r 2i~0; (1.72) 

£ sl a4-Ar M ?+7’ al =-0, (1.73) 

> 

— over the section after separation of the first stage 

•5 4 =yi s2 o 5 + ^“+A r s! H7' M =0; (1.74) 

5 5 =^ 32 a 2 +£ 32 a + Af 3a 3+r3 2 =0. (1/.75) 


Here 
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Let us present an example of the selection of the optimal con- /43 
trol a and $ at some fixed instant of time t, where we will assume 
that at this instant no staging occurs. The selection of the opti- 
mal control at the instant of staging will be considered a little 
later. 

As we can see Eqs . (1.70) and (1.71) are identical with 
respect to variables a and 3 and differ only in the coefficients 
appearing by these variables. These coefficients at any instant 
of time (except for the instant of staging) are uniquely defined. 

As a result, to condense the notation, we replace Eqs. (1.70) and 
(1.71) with a single formula of the form 


H=Aa i +Ea + Fat-\-Wt-T. \ 


( 1 . 76 ) 



The coefficients A, E, F, and F of these formulas are equal to the 
following : 

— over the section before separation of the first stage 

A=A n ; E=E U \ E=E n ; N=N u ; T=T n ; 

— over the section after separation of the first stage 

A = A 12 \ E~E xz \ E = E n \ N = N l2 \ T~T n . 

It is required to find in Eq. (1.76) a and 3 satisfying the 
inequalities 

a 2mln ^ ® ^ a 2ma*’ ) 1 

' Pmm«P<P,„«. I ! (1.77) 


where the function H would take on its largest value. 

Let us examine the case when the constraints (1.49) and (1.50) 
have been satisfied. We will show that the function H. takes on 
its largest value only along the boundaries 3 = 3 m in and 3 = 3 ma x' 

Let us set up the determinant: 

dm dm 

da 2 dad$ 

om dm 

dad$ d$i 


This shows that the function H can have within the specified 
region [domain] (1.77) only a saddle-point . 

Since the mass flow rate of 3 appears linearly in the function 
of H, this means that along the boundaries a = 012 m i n and 
a = a>2 max ? the function H takes on its largest value either when 
3 * 3 m in or when 3 55 3maxj that is, the maximum of the function H 
can be realized only at the apices of the rectangle (1.77). Thus, 
the maximum of the function H can be realized only along the 
boundaries 3 - 3 m in and $ “ Smax- 

We will present points at which the extremum of the function /44 
(1.76) can be realized as the LV moves within the domains (1.49) 
and (1.50), in final form: 


2 A F 
F 0 


— — < 0 . , 


(1.78) 
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1) 

a 2ni l n 



2) 

a 2max 

n 

• m I n ’ 


3) 

a 2in 1 11 

r'max 1 


4) 

a 2max 

a 

.max* 


5) 

= 


' ft 

2A •‘ n1 "’ 

6) 


E 4 

■ £?max ft 

a* — 


a nnax ’ 

2A 


(1.79) 


To condense the notation for defining the optimal control as 
the LV moves along the boundaries (1.49) and (1.50), let us 
introduce some polynomial 


y = £ ft a*+ B t a+ fl,=0. 


( 1 . 80 ) 


The coefficients Bq> Bp, B 2 , and B 3 of the polynomial are 
t&ken as equal to the following: 

over the section before separation of the first stage: 

a) for motion along the boundary £ 4=0 

B 0 = An x ‘, B 1 ~ B 2 — N ai : ^3—^211 

b) for motion along the boundary £5 = 0 

^0 ^31' B 2 — N 3 i ; B 3 ~T 91 ; 

along the section after separation of the first stage 

a) for motion along the boundary £4 « 0 

! 

5 0 = vA 22 ; B y = E^ B Z =N M ; 

b) for motion along the boundary - 0 

Bq — A 32 ,', B 1 ~E 32 \ B<z — N B 3 —T 3 3 . 

Polynomial (1.80) is a generalization of the polynomials 
( 1 . 72 ), (1.73), (1.75), and ( 1 . 75 ) we have already introduced. 

By determining & from this polynomial and inserting it into 
Eq. (1.76), we will have 


// = *4<* 3 -f B 5 a 2 -f B s a.+ B„ ! (1. 8 l> 
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where 


B t - B,= A-~ '-(FBt+NB,,); 

A IB* 


B^ — E — ~{F Bs+NBJ, B 7 ~T 
»■> 


B 2 


Then we can quite readily find the points at which the func- /45 
tion H can have an extremum for motion of the LV along the bounds 
(1.49) and (1.50)* They will be: 


— B 5 — \/~ B\ • — 35^6 ^0 a 3 + ^l a 3 + ^3 

1 ) a 3 — - = 


3 B 4 


Bn 


B\ + t/* B\ — 4 B 0 C^Praln + ^3) 

2 ) o,= 

- B, + V B;~4BoWn.i« + r 

3 ) a 5 ~ 


2 B„ 


4) u e 


— B\ + J/ Bj — 4B 0 (B 2 ^, nax +* B 3 ) 

5) a 7 — „ Pma*’ 


®) a !imln 
1 ) ® 2 max 


1^2 — 


ft — _ 

^3 


2 B 0 

jg 0 o.| n , j „ + B]Ct2 mt ,i 4- B 3 
fl 2 - 

fi-l g 2max + ffl a 2max + B3 

bT” 


Bj .{- Y B l~ 4jB f’ (^Pnmx + fi 3) , } 

— | 


( 1 . 82 ) 


The values of a^, a 2 , cto, a/j , ac, ag.* and ay, and also 81 , 8 2 s 
and 83 obtained by Eqs. (1.79) and (1.82; must satsify the cor- 
responding conditions ( 1 . 77 ). If for each point these conditions 
are satisfied for both a and 8 , and, in addition, the conditions 
(1.49) and ( 1 . 50 ) are satisfied, the resulting point is left for 
further investigation; but if these conditions are not satisfied 
or if one of the values of a or 6 are imaginary, it is discarded. 


It should be noted that at any instant of time the number of 
points satisfying these conditions is smaller than is presented 
here. Actually, many of these points preclude each other. Some 
of them are imaginary. However, we cannot in advance say exactly 
to which . point at a given instant of time the optimal control cdr^ 
responds:,. ..Therefore, we must check all these points and leave 
those which satisfy our requirements. 
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In the following we will not qualify this in every instance, 
and if we consider a particular point, we will assume that the 
conditions specified for it have been met. 

By calculating the values of the function (1.76) for points 
satisfying all the above-listed requirements, let us select the 
point at which the function H takes on its exact upper bound. The /46 
control a and 8 thus obtained is the optimal control at any fixed 
instant of time t (except for the instant of staging). 


Determining Optimal Control at Instant of Staging 

As already indicated earlier, at the instant of staging $min* 
Bmaxj 5, P S p, c Xq> a min> and a max can suffer discontinuities, as 

can other LV parameters; the Lagrangian also suffer discontinuities. 

We are required, speaking figuratively, to "dovetail" two 
different sections: the section prior to staging and the section 

after staging. Condition (1.59) will be the criterion of this 
"dovetailing," and this condition must be realized for an exact 
upper bound of function H + in terms of control. 

Prom conditions (1.62), we have: 

( i== 1 * 2 » 3 , 4 ), 

'*5 + =*r+ v * (1-83) 


that is, at the instant of staging only the multiplier Ac suffers 
a discontinuity, and the rest are continuous. 

The function H + at the instant of staging is of the form 
(1.71). All the coefficients of this function, with the exception 
of Nq2? are known at this instant, since the characteristics of 
the LV after staging are specified ahead of time, and the multi- 
pliers Aj_ (i = 1, 2, 3, 4) at the instant of staging are con- 
tinuous. As for the coefficient Nq2> because the multiplier 
has a discontinuity, we as yet do not know this coefficient. 

Let us note that the formula functions (1.79) and (1.82) 
serving for the determination of the optimal control within the 
sections separated by surface (1.54) are valid also at the instant 
of staging, and in those functions where the coefficient N does 
not appear, this control Is uniquely defined both at the left and 
at the right, and it requires only that into these functions the 
LV parameters corresponding to the given stage be substituted. 
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As we can see from Eq. (1.79), the coefficient N does not 
appear in the formulas for the optimal control, and in the func- 
tions ( 1 . 82 ) it appears only in formula 1). 

Let us define the optimal control and thus the coefficient 
N]_2 for the case when staging occurs within the domains (1.49) and 
( 1 . 50 ). 

Let us write the condition (1.59) at the instant of staging 
in expanded form: 


H~ — A n a~ 2 J - F n a -j- F u a '-i -}- A 7 n p -{- T u — 
— A l 2 a + ' 2 -j~ £’ ls cr -f ^i2 a+ ? + + N ja ? + — 


Since H at the instant of staging is known, condition (1.84) /47 

can be represented as: 


■ 4 1 2 a+2 + ^r 3 a "M“ T^ — l2 a + + A^ 12 ) [i + =0. K±-ooj 

We note that only one of the limiting values 3 m i n or 3 ma x 
can be the optimal value 3 + , just as 3“ (with the only difference 
that their numerical values [ 8 m in and 3max^ correspond to the 
second stage and can differ from the numerical values of 3 m in and 
Smax at the first stage). 

Let us denote the part of the "function (I. 85 ) independent of 
3 + by C, that is 


C...A a <x'" j- E^+Tn-H-. '■ (1.86) 

Let us substitute here the values a + corresponding to points 
1, 2 and 5 from formulas (1.79). 

We note that the value of 3 + corresponding to point 5 is 
substituted into formula (1.86) only if it satisfies condition 
(1.77), and the coefficient A 12 <0. It is not considered in the 
remaining cases. 

Since the points 1, 2 and 5 correspond to 3 = 8min> this 

means that the following must obtain at these points: 


f’i2 a ' h -f w 13 <o, : 

or, because 3 m ih > 0, C must be greater than zero. 


(1.87) 
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Let us select from points 1, 2, and 5 those for which C > 0, 
and let us reject the remaining. 

Let us substitute the remaining points into the formula 

; ( 1 . 88 ) 

and, let us dwell on the points for which Ci is at a maximum. The 
coefficient N12 corresponding to the points thus selected is 
defined by the formula 


yy H ' ~ ^- Q+2 ~~ £l20t+ ~~ /r '> a *P + ~ 7 ‘ 2 . ' (1.89) 

12 ?'* I 

Let us substitute into Eq. (1.86) the values a + corresponding 
to points 3, 4, and 6 of formulas (1.79), where the latter point 
is substituted into this formula only if the conditions (1.77) 
and A12 < 0 are satisfied for it. Let us select from these points 
those for which C < 0, and let us reject the rest. 

We substitute the remaining points into Eq . (1.88) and we •' 
dwell on the point at which has a' maximum 1 .” The coefficient N]_2 
corresponding to this point is defined by formula (1.89). 

By carrying out this study, we find at each of the bounds 
3 - 3 m in and 3 = Bmax one point each (the case when there is no 
point of this kind along any of the bounds is possible), at which 
the exact upper bound of function H + can be realized. 

Now we must select from these two points the one for which, /48 
at the coefficient N12 found for it the function H+ realizes its 
exact upper bound at this point. 

Obviously, if we take as the optimal control the point for 
which the coefficient N12 is the smallest, we obtain the result 
that at any other point for the coefficient Ni2 thus determined the 
function is smaller than at the point we selected. 

Thus defining the optimal control a + and B + and the coeffi- 
cient Ni2, we can also find the multiplier X^ + : 




(1.90) 


Now let us show that the control thus selected always exists. 

If A]_2 < 0 and a + a 2 + (see Eq . (1.79)) satisfies the condi- 

tion (1.77) > by substituting and ot2 + into Eq . (1.86) we can see 


43 



that C(a-]_ + ) > C(a2 + ). Hence it follows that if C(aq + ) < 0 (the 
case when the optimal control is absent along the bound 3+ = 8 m q n , 
then C(a 2 + ) < 0 and this means that the optimal control exists 
along the bound $ + = Bmax. If C(a 2 + ) > 0 (the case when there is 
not optimal control along the bound B + = $max) * then C(aq + ) > 0 
and this means that the optimal control exists along the bound 

~ ^min* 

If A]_2 > 0 (or A12 < 0 > but aq + and a2 + do* not satisfy condi- 
tion (1.77)), obviously the function H + realizes its maximum along 
the boundaries a + - a2 m ± n and a + - a 2 max • Substituting into 
Eq. (1.86) the values of a+ , we uniquely define the sign of C, 
and hence thus the values of B + corresponding to the bounds. 

Therefore, in this case there is an optimal control. 

Now let us assume that A12 < 0, and that one of the values of 
aq + or a2 + (for example, aq + ) does not satisfy condition (1.77). 
Obviously, if a]_ + > 012 max , then we must select a2 mix ; butt if 
a l + < a 2 min 3 we must take 012 min • Let us denote! this control by 
the symbol ag. Since a2 + satisfies condition (1.88) and aq + does 
not, then clearly ag is in this segment [aq + ,ag + ]. On the other 
hand, the maximum of the function C (see Eq . (1.86)) is realized 
when as ““Eqp/^Aqp* By comparing the formulas for aq and a 2 + 

(see Eq. (1.79)) with the formula for a8 , we see that the value of 
a8 lies outside the segment [aq + ,a 2 + ]. This means that in this 
segment the function C either decreases or increases with respect 
to a. We earlier showed that C(aq + ) > C(agt). Now we can write 
that C(aq + ) > C(an) > C(a2 + ) . Further considerations are analogous 
to the case when A12 < 0 and aq + , a 2 + satisfy condition (1.77). 

Thus, at the instant of staging the optimal control always 
exists and does so uniquely (an exception may be found in several /49 
singular cases, which we do not consider here). 

Now let us define the optimal control and thus the coeffi- 
cient Nqp f°r the case when staging occurs along the bounds (1.49) 
and( 1 .50} . 

As already indicated, in the formulas D'of the function (1.82) 
there appears the unknown coefficient N]_2 . So aq and thus 81 
cannot be determined directly from these formulas. 

We can see that at the instant of staging the coefficient 
B4+ and do not depend on the coefficient N12 , while Bg + and 

By + are functions of Nq2 • us represent Bg , Ni2, and By + as 

functions of the angle of attack 013+, and the coefficients 
Bi| + and B5 + (this can be done by using formula 1) of functions 
(1.82) and the coefficients earlier obtained for the formula 
(1.8D). 
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( 1 . 91 ) 




3fi 4 + a 3 +2 + B+ 


2a ; t 


N 12 — 


3 (B+a 3 ' 2 + 2A[&1' + B <l) B 2 , 


Bt~T l% - 


2B+at | 

1 2/1^+^^)^' ! 


2B„ w 3 


(1.92) 

(1.93) 


Substituting the resulting expressions for B^ + and By+ into 
Eq. (1.81) and considering condition (1.59) * we have 


where 


B. 


4 -J- /^ 0 a:! -+ /*,.»«» I 
fif. fi w --=2(7V #«' ’ 






BtBt 


(1.94) 


In this expression coefficients Bg , Bg , B^g> and B^n are 
uniquely defined and the problem thus reduces to 'tfie solution of 
this equation in a^ + . To solve it, we can use, for Instance, the 
method of tangents or chords [ 1 8 ] . In solving this equations we 
must find only the values of 013 + that satisfy conditions ( 1 . 77 ). 

By finding and then the corresponding to it, and by 

verifying that conditions ( 1 . 77 ) have been satisfied for the 
values of cto + and Bi + thus obtained, let us dwell on the values for 
which conditions (1.77) are satisfied. If there are no such 
values, then let us move on to considering the next points in 
formulas ( 1 . 82 ) . 

Let us determine for the values a^ + and 8 q + satisfying condi- 
tions ( 1 . 77 ) the coefficient N ]_2 corresponding to it based on 
Eq. (1.92). Then, let us substitute the values of a 3 + into /50 

(1.86) and let us define the sign of the function C. If C < 0, 
ct 3 + and Bi + are left for further investigation, but if C >0, they 
are discarded. But if after all this analysis it turns out that 
there are several solutions to Eq. (1.94) for which all the above- 
listed conditions are satisfied, then let us dwell on the solution 
for which the coefficient N 12 Is at a minimum. 

The value of a + and B + for the remaining points of formulas 
(1.82) at the instant of staging are determined directly. By 
verifying condition ( 1 . 77 ) for these points, we discard the point 
for which these conditions are not satisfied. For the remaining 
points we find the sign of function C (see Eq. (186)). For points 
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2 and 3 we must have C > 0, and for points 4, 5, 6 , and 7 there 
must be C < 0. By verifying this analysis, we discard the points 
for which the sign of function C does not correspond to the re- 
quired sign. 

Let us find for the remaining points the coefficient N ]_2 
corresponding to them, using the formula 

; v gjT «q +a + 'M * 2 + K a+ +T\2 -h~) ( 1 . 95 ) 

12 -fi n + a+2-f B+ 

Knowing Np 2 * we can find the coefficient By + of Eq . (1.8l). 
Obviously, if we take as the criterion of the optimal control the 
control for which the coefficient By + is the smallest, the 
function ( 1 . 8 l) will take on the largest value with respect (-to 
the control, that is, the control will be optimal. By thus finding 
the optimal control and the coefficient Np 2 corresponding to it, 
we can find the multiplier X^ + by formula ( 1 . 90 ). 


Solution of the Problem 


In the preceding sections we outlined a method for determining 
the optimal control at a fixed instant of time x. By thus deter- 
mining the optimal control, we can now find the multipliers ^ 3 , 
ij^, and appearing in this system of equations (1.64). 

They are defined from the conditions 

2,4 lft a + £ lk + « 0 , ( 2 a - a m , n - 0 + 

c a y Sqr 2 ' i 

H" 'W ^yiraln n in max J “ T i j (l 96) 

4- *-l» 4 (2A 2i a "f £’as)+ ^6 (^ 3 * a + £ 3ft); J I 

F Jk a 4 -Njjt — (2? — ? ra i n — «) 4 'KN 


Here k = 1 and k * 2 after staging. 


As we can easily see, the two equations contain five unknowns. /51 
Since there are two control functions, motion is possible only 
either within the domains (1.46) - ( 1 . 50 ) or else along one or 
two of the bounds of these domains. Hence it follows that the 
maximum of the multiplier (j = 1, 2, 3, 4, and 5) simultaneously 
nonzero cannot exceed two. 


4l 0 , 
j- 0 , 
^3 ¥ 0, 
^4 ¥ 0 , 
^5^0, 


if for the selected optimal control 

it n 11 u n u 


it n if 

it n it 

H it ti 


Cl = 0 ; 

C 2 = 0; 
C 3 - 0; 
Cii = 0 ; 
C 5 - 0 . 
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The multipliers are equal to zero in all the remaining cases. 

Next* let us integrate the systems ( 1 • '4 5 ) and (1.64) with the 
optimal values of a and 3 thus obtained, and also i/u . In the 
next stage of integration the entire process for determining the 
optimal controls a and 3 , and also the values of is repeated 
over again, and so on. 

The conditions for the discontinuity of the multipliers X-^ 

(i = 1 , 2 , 3, 4, 5 ) at the instant of insertion are as follows: 


— at the bound ^4 


0 

X+=a r -f-v^; 


‘' 3 ' 7 ^ ’ 

X;= X- + V J2.!2.; 

4 4 Ok 2 




1 

1 

1 

I 


(1.97) 


— at the bound ^5 


1 1 ' 

K i= v; 


v.: 


x* — v , d li*L- 

4 4 1 Oh 

A+ = X“ 


(1.98) 


Now let us formulate the general solution to the problem. 

Since the function H is homogeneous, the initial value of one 
of the multipliers Aj_ can be selected arbitrarily, for example, 

^•50 = ^ (1*99) 

One of the multipliers X i9 for example ,> X 40 , can be determined 
from the condition: 


x 


40 


X10/1 -t h‘J\f 2 - 1 - hip/a + X 5[) /s 
f 4 


(1.100) 


When the systems (1.45) and (1,64) are integrated, we have /52 
ten arbitrary constants, and the time of completion of insertion 
tf is also unknown. 
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To determine the 11 unknowns, we have 11 relations ( 1 . 52 ), 
( 1 . 53 ), ( 1 . 99 ) and (. 100 ). Integrating the. systems ( 1 . 45 ) and 
( 1 . 64 ) under the initial conditions (1.52), ( 1 . 99 ), and " ( 1: 100 ) , 
and using, for example, Newton T s method [ 18 ] to determine A^q, 

^20i and A3Q, we satisfy the conditions ( 1 . 53 ). 

In solving the problem with constraints as on the first-order 
phase coordinates, the initial value of one of the multipliers Af, 
for example, Aqo, is selected so that condition (1.59) is satis- 
fied at the instant of insertion at the bound £4 = 0, and then by 
selecting the number v, let us satisfy the conditions (1.59) at 
the instant of insertion at bound 65-0 and the number Vf 
(equally with A 2 o and A30) is determined on the condition that 
the conditions (1.53) are satisfied. 

First, let us dwell on the solution to the boundary value 
problem without referring to the constraints on the phase •<•.! r ' 
coordinates . 

Suppose that in the first approximation of the initial multi- 
pliers Aqoj A20 3 and A^q we assumed A 1 -j_ 0 , A 1 20 and A *30, and that 
the resulting mismatches of the boundary conditions corresponding 
to this approximation are as follows: 

W — V — 

A& 0f, 

| 

a! —L — Zf , \ 


where V, 0 , and L are the valuesobT V, 0 , and L obtained by in- 
tegrating the systems(l. 45 ) and ( 1 . 64 ) at the instant h = hf. 

Varying A’^O , A* 20 3 and^A'go by a sufficiently small quantity, 
let us find at the instant h - hf the partial derivatives 


ov ov OV <90 dti 
— r> — * » ’ > 1 . . 1 

^30 ^10 ^ a 20 



di 01 of- 


Knowing the partial derivatives, we can compute the correction 
AA’qOs AA T 20 > and'AA r 30 ,to the Individual values A’qo, A T 2 0 3 and 
A' 30, for example, by Newton r s method, thereby solving the -system 
of hfchree^algebraic.linhomogeneous equations 


dV w , 0V 

-rV 10 i r 

ffk\« oi 


A^20“T 


10 

dXi 0 

0L 


1^0+ 


20 

dfl 


aa;- 1 - 


OV 

^30 

d0 


AA 30 -t-sA^ — 0;, 


10 1 - » ■ • . AX; o +eA 6 - 0 ; j. 

OKr? 01,- ' 


; A a jq 1 

dx; 0 oi 


20 

dL ,, , dfl 
A i" 


m i 

AA' -}- eA L~ 0, 


20 


0 \ 


30 


) 


( 1 . 101 ,) 
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The process is repeated until the required accuracy is reached. /53 

The quantity e (0 < e <_ 1) is selected on the condition that 
the solution of the boundary value problem converges, where the 
closer we come to the assigned boundary conditions, the closer 
£ is taken to 1 . The speed with which the optimal solution of this 
problem 1 is determined depends on the proper choice of e. 

If the first approximation is sufficiently good, we can assume 
e = 1 . But if it is far from the required value, the value of 
£ must be taken closer to zero, by increasing its value from one 
approximation to the other. 

If the problem is solved with reference to constraints on the 
phase coordinate, the process of solving the boundary value prob- 
lem is somewhat more complicated, since in this case, in addition 
to everything else, the conditions of discontinuity at the docking 
moment must be satisfied. In turn, the new multipliers v and vq 
will appear, which also must be determined from the solution to 
the boundary value problem. 

To solve this problem, it is useful first to solve the prob- 
lemtwithout reference to the constraints on the first-order phase 
coordinates. Then, by assuming to the first approximations the 
resulting values of X]_q, X203 and X30S and revising Xqg, X2g, X3Q, 
v, and v^, we satisfy all the conditions of discontinuity at the 
moment of docking and the boundary conditions at the right 
endpoint of integration. 

Let us consider the solution to this problem more closely. 

Suppose that we adopt as the first approximation of the 
initial multipliers Xqg, *20* and X30, the quantities X’qg, X^Oj 
and X'3o, and suppose that the function H+ is not equal to zero 
at the moment of docking. By varying X*qo by an sufficient ly* smalls 
quantity, at the moment of docking we find the partial derivative 

9 H + 

3 x • The multipliers XX 2Q and: X ! ^g-.are 'not yet Varied. 

By knowing the partial derivative, we can compute the cor- 
rection X*qo for the initial value X'qg based on the formula 


’ 10 "" ''dH+Zdl^ 

The process is repeated :until the required accuracy in H + is 
attained. The value of v is revised similarly, only it is carried 
out now on the condition that the< discontinuity is satisfied at 
the moment of second docking, namely, by assigning the approximate 
number vq at the instant of first docking, we determine H+ at the 
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moment of the second docking, therefore varying v °\ by a suffi- 
ciently small quantity, let us determine ! ;the partial derivative:. 
3H+ 

3 ^ at the moment of second docking. Then the correction to the 
initial value vo will be determined by the formula 

H + ' 

Av tt — — c a dli +/dv Q ‘ | 


The process is repeated until the required accuracy in H + is 754 
attained at the moment of second docking. 

The value of Vi is revised on the condition that, when h = hp, 
condition Y = Vf Is satisfied, based on the formula 


Av 10 —- 


AV 

3 oV/dv l0 ‘ ■ 


The trajectory thus selected satisfies all the conditions of 
discontinuity at the moment of docking and the two boundary condi- 
tions ( 1 . 53 ) at the right endpoint of integration. 

Then, by varying A^O and A'oq by a small enough quantity and 
by solving the analogous problem In the selection of Aqo> and 
vq , let us determine the partial derivatives 

do d6 OL OL 

r j , > , * , * , 

dX 20 dX 3Q dX 20 c/a . 30 


Knowing the partial derivatives we can compute the corrections 
AX 1 20 an( * AX T 30 to the initial values A^o and X r ^o 3 by solving the 
following system of two equations: 


do 

dXoQ 

OL 

dX 2 o 


A X 20 "T A X 30 + £ 4 — Oi 

.^30 

AX'o + A X 30 + = °- 

^‘30 


I 

I 

I 


The process is repeated until the required accuracy in satis- 
fying the boundary conditions ( 1 . 53 ) is attained). 

The values of the numbers eq, £2> e 3 and £4 ( 0 << eq <_ 1, i = 

= 1 , 2, 3, 4), 'Just as e, are selected on the condition that the 
process of solving the problem converges. 
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Particular Cases of the Overall Problem of Optimizing; the Tra ( jec4 
to ry of LV Insertion into AES Orbit 

First let us examine the case when an LPRE is used as the 
power plant in both the first and second stages. We know that the 
optimal trajectories of an LV with LPRE areocharacterized by its 
rapid departure from the dense atmospheric layers, which enables 
us to remove the constraints (1.49) and (1.50). Thus, Eqs; (1.56) 
and (1.57) will be of the form 


x 

OH , 

Aj 


dxi 1 

^ OXi 

OH 

f'l 

dSi 1 


~r~ — Ti 
0u 

Ou ' 

1 “ Ou 


(/=!, 2, 3, 4, 5); 


? 3 -r («*= .*)■ ; 

ou 


( 1 . 102 ) 

(1.103) 


Condition (1.59) in this case (with the exception of the mo- /55 
ment of staging) will be automatically satisfied, and conditions 
(1.60) and (1.6l) will not be used. 

Eqs. (1.64) can be written thusly 


•; \ d/i x .'d/3 a i • 


2 1 $8 2 ^ 3 


OV 

0f2 

ae 


OV 
Of a 
d(t 


ov 
Of A 

- dO 






a/L -h4— v ^+ t. 


dh 
Of 1 


dh 

df<z 


Of* 

dh 

dH; 


db 

dh 


i _ 1 <*J\ x J l I & 

^6 — h A 2 , I T, it 1 

dm Orn 0m 


(1.104) 


Of 1 Of 2 Of ‘A Of 4 
The partial derivatives ( ,v ' ov ’ ov 1 ov ’ 


Oli Of I Of 2 
ov ' do * as 



; d/4 Of i Of 2 
\ d(l ’ dh ' dh ’ 


d/3 , d£ 3 Of 1 

dh ’ dh ' dm 


- are determined from the 

dm dm 


same formulas as in Eqs. (1.64). The engine thrust in both the 
first and second stages will be determined by Eq. (1.67), and 
its partial derivatives will be of the form 


OP 

ov 


0, 


OP 


~go e»S a 


dQ 


dh ^ “ o fl 

But the formulas (1.65) and (1.66) are not used. 


(1.105) 


The coefficient Ail, E^, F^, Fn, Nn and Tn of Eq. (1.70) 
will be of the form of the corresponding coefficients A 12 , Eqg , 
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Ft_ 2 » and T 12 of Eq. (1,71). Their only difference will he 

that in the computation of coefficients with subscripts " 11 " we 
must take all the characteristics of the LV and the power plant 
corresponding to the first stage, while when computing the co- 
efficients with subscripts ” 12 " we must take the corresponding 
characteristics of the LV and the power plant of the second stage. 

Because the constraints (1.49) and (1.50) are removed, the 
case of considering the motion of an LV along the bounds £4 and 
£5 disappears and, naturally, the formulas ( 1 . 72 ), (1.73), (1.74), 
(1.75), ( 1 . 80 ), ( 1 . 81 ), ( 1 . 82 ), (1.91), (1.92), (1.93), (1.94), 

( 1 . 95 ), ( 1 . 97 ), and ( 1 . 98 ), just as the solution of the boundary 
value problem with reference to the constraints on the first-order 
phase coordinates, will become unnecessary. 

In the equations (I. 96 ), considering the foregoing, we must 
adopt 4 55 $5 * 0. All the remaining formulas, which are not 
mentioned here, are used In the form in which they are represented 
in the method given earlier. 

In conclusion, let us dwell on the insertion of an LV into 756 

an AES orbit, beginning with the end of the atmospheric section 

of the flight when LPRE are used in the second and third stages, 

and adopting as the specific thrust the specific thrust of the 
LPRE in vacuo. 

In. this case, system (1.45) becomes: 


y_2^_JL s in9— / x ; 

m r- 

a- . — e_ cos 

mV r*V r 

iL/cosB-/,; 

r 

fi V sin 0 = / 4 ; 
w— — P = / 6 * 

Naturally, In this formulation of the problem we must con- 
sider only the constraints (1.46) and (1.47), while we must 
discard the constraints (1.48), (1.49) and (l.$0). 

The system of conjugate equations ( 1 . 56 ) in this case will be 
of the form: 


— fv 


(1.106) 


A; 


OH 


dxi 


(i— 1, 2, 3, 4, 5). 


(1.107) 
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Let us write it in expanded form: 


A, - 


* V m\ 


? a 


u 0 COS 

cos 6 ■ 


r } r 


cos 0 — X 4 sin 6 ; 


iV 2 , 21/2 

^a — ^i — cos 0 — ^3 ( — — ) sin 0 + — V sin 0 — \V cos 0; 

/• 2 “ \/-2K r J r 

>•3=0; 

^=- i i-^-sin6-i s ^Jt.-ilJcos() + i 3 -a KcosO; 

i. = 4^ (*,+*,£)• 

in 2 V f 


( 1 . 108 ) 


The function H both before separation as well as after sepa- 
ration will be of the form 


H = Fa$-\- N$-\-T> (1.109) 


where 


mV ’ 

* = *.- 2I»_X 6 ; 

7 " - ** 7 T sin 9 + ** (7- ' ^7 ) cos 6 

+ >'3 ~ V cos 0 4- X d l/ sin 0 . 
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Since — y ■ > 0, the sign of coefficient P is uniquely deter- 
mined by the sign of multiplier X 2 • Since 3^0, it is clear that 
the sign of F3 will be characterized by the sign of X 2 • And since 
a appears linearly in the function H, the optimal control in terms 
of the angle of attack a can take on!; only limiting values, where 
if X 2 s °* we must taken a opt = a mi n; if A 2 > a opt = ct max . 

Similarly, we can show that only the limiting values defined 
as follows can be optimal controls in terms of the mass flow rate 
of fuel 3: 3opt ~ 3mina if Fot + N < 0j 3opt ~ ^maxa if 
Fa + N > 0. 

Thus, the optimal control in any fixed moment of time x (ex- 
cept for the staging moment) is uniquely defined. 
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Now let us dwell on selecting the optimal control at the 
instant of staging. In the case examined here, the problem of 
determining the optimal control at the instant of staging is solved 
quite readily. 

Actually, at the instant of staging only the multiplier X 5 
suffers a discontinuity (see Eq. (I.83)). This multiplier appears 
in the function H in terms of multiplier N. Since the angle of 
attack a is entirely defined by the sign of multiplier X 2 (and 
this multiplier at the instant of staging is continuous), it is 
clear that the angle of attack a + at the instant of staging is 
determinedcuniquely based on the method indicated above. 

By computing the coefficient T + (we know that T + - T~ , since 
the rocket parameters do not appear in the multiplier T), we find 
its sign. T + < 0, then because H + = 0, it follows that 
(F + a + + N + ) 3 + > or expressed another way: 

F + a + + N + > 0. 

And this means that we must take B + = 3 max as the optimal 
control at the instant of staging. 

But if T + > 0, we must take 3 + = 3min & s the optimal control. 

The coefficients N + and X^ + corresponding to the control thus 
selected are determined by the formulas 


N + — 


7 + 4 - F i a+ 




J ° go . __ 7V + . 
m~- 


1 , 3 . Method of Inserting Satellite into Orbit Without Final 758 

Combu s t i on "of Fuel 


Let us examine the case when a satellite is inserted into 
orbit with a continuous powered trajectory section of the launch 
vehicle. Obviously, the shape of the orbit and its distance at 
different points from the Earth's surface will depend on the 
selection of the elements of the end of the powered section, and 
also on the adopted law of launch vehicle motion. For many kinds 
of; satellites it is important that their orbital perigee lie as 
high as possible. Accordingly, we will seek the maximum altitude 
of the end of the powered section for an assigned final Velocity 
and assigned direction of this velocity. 
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We will not dwe-li: on methods of calculating the powered sec- 
tion itself, assuming it to be known to the reader from already 
published studies (for example, [5]). 

Let us determine the optimal law of launch vehicle motion 
along the powered section, that is, let us find the optimal pro- 
gram of variation in its pitch angle <j>p r , without referring to the 
variable fuel consumption. The methods examined below arena 
further elaboration of the study [46], They permit the quite the 
quite simple selection of 4>p r virtually without loss of accuracy. 
The optimal 4>pr found can be used in an exact calculation of the 
powered section for a specific launch vehicle. 


Optimal Program for Angle of Pitch with Constant Fuel Consumption 
over the Powered Section 


We will start from the system of differential equations of 
motion that take into account the curvature of the Earth and the 
variation in gravity with* altitude [46]. For the case of constant 
fuel consumption this enables us to obtain a simple analytic 
method of calculation: 


where 


(l U .1 

^ p cos — vtr; 

dt 


dw 

dt 

dx 

dt 


— p sin cp -f ‘2v 2 y — 


u; 


dy 

w, 

dt 


( 1 . 110 ) 


v 2 = 



go is the acceleration of the force of gravity at zero ■ altitude; 

R is the Earth T s radius; p is launch vehicle acceleration; and 
4) is the angle of pitch. 

The system of equations (1.110) can be readily obtained from /59 
the ordinary equations of rocxet plane motion written in an ab- 
solute rectangular system of coordinates with the origin at the 
launch point, the Oy axis directed upward along the vertical, and 
the Ox axis directed toward the powered section of i the trajectory. 

By expanding in power series the expressions for the pro- 
jections of gravitational accelerations 
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gy = go 


R*x 

[ X 2 + Q,+ R) 2] 3 / 2 ’ 

R? (£ + R) 

[ x * 4 - (« + ^) 2 ] 3/2 


I 


and by limiting ourselves to the first terms in the expansion, we 
get 


„ X i 

gx — 1 

gy — go'^'^go ; 


In the particular case when we are considering a two-stage rocket, 
we will take as the initial point the end of the powered section 
of the first stage, designating it by the subscript ''l". 

Let us set up the variational problem of finding the function 
4 >(t) for assigned boundary conditions that give a maximum to the 
altitude 



( 1 . 11 ) 


where the subscript n f" refers to the elements of the end of the 
powered section. 

It should be noted that this problem is equivalent to the 
problem of determining the function <J>(t) that gives a maximum for 
the velocity at assigned altitude [ 46 ], From the standpoint of the 
calculus of variations, this is a problem of the conditional ex- 
tremum with differential relations. The equations of motion 
(1.110) of the launch vetiicle are the relations. 

Let us set up the Lagrangian function 


H -j- Ag'jjg -j- A 4 '|» 4 , 

where A^, Ag , A3, and A4 are certain differentiable functions with 
respect to time t: 
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4. — — A-v^x— jO cos cp; 

' 1 at. 

= 2vhj - p sin <p -f g 0 ; | 


The functional whose maximum we seek is of the form 


=»A f + ftfrf 


( 1 . 112 ) 


Let us write the first variation for I: 


S' — SA:f-f~ | 3 Hdt , 


(1.113) 


S// i- "77“ : 

1 R i 

3// - H u %u -f H w ^.w \~ H \jly f - H x %x 4 H „ bu 4 
-r 4- H 4 Hjhx 4 H »8(p. | 

Canceling out the derivatives 6u, 6w, Sy, and 5x by integrating 
the corresponding expressions by parts, we get 

8 l = ?,y f +^- 8*f+ |>-,M' f + [W' f ■ + ; 


+ W</|' £ + f 

6 

i 

• ■■{'2v 2 l 2J r ). i )<iy 4 pO-i sin ? ■ - cos <?) 8<p] (it. \ 


The extremal value of the functional I corresponds to the 
zero-equality of its first variation 61, 

By setting equal to zero the multipliers 'appearing under the 
sign of the integral in front of the variations 6u, 6w, 6x, and Sy, 
we get a system of differential equations for the quantities 
A 2 > and ^4 and the final relation for <f> : 
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(1.115) 


— 0; I | 

X 2 -|-X 4 =0; [ | 

X 8 — v 2 a 2 =0; j 

X 4 -{-2v 2 X 2 = 0; 

X x sin — X. 2 cos cp=0. i 

From the last equation we have 

I 

f x s i 

The solution to the system of equations 

form: 

X t = C x cos vt C a sin vt; 

>^=C 3 sh /Sv/-fC 4 chK2^; 

X 3 =C 1 v sin W — C a v cos v*; 

X 4 = -C 3 V"2vch /2vif-C 4 /2vsh]/2v/. 

Eqs. (1.115) and formula (1.116) are analogous to the cor- 
responding expressions in [46]. Let us select the following 
boundary conditions. At the end of the powered section we will fix 
the velocity vector (uf , Wf) , and at the initial point — the 
modulus of velocity, that is, at the end of the powered section of 
the first stage the following relation must be satisfied 


(1.117) 


(1.116) 

(1.115) is of the /6 1 


vf uf | wf const, 


(1.118) 


which gives for the variations 6u^ and Sw]_ the relation 


i (1.119) 

«i i 

We will also fix the coordinates x and y of the initial point. 

This selection of the boundary conditions simplifies the solu- 
tion of the problem. It is quite justified, since assigning a 
constant velocity modulus at the initial point instead of the 
functional dependence of velocity on altitude and angle of inclina- 
tion of the orbit Is based on compensating for the effect of change 
in velocity at the initial point on the insertion altitude h f by 
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the influence of the corresonding change in altitude (in angle) 
at this same point. 

It is convenient to take as the initial point .the < end off 
the powered section of the first stage, since in the first stage 
of flight occurring in the dense atmospheric layers, 4> pr is 
selected on the condition that the angles of attack are near-zero, 
and not on the condition that the maximum altitude of the insertion 
point is obtained. 

Let us examine the case of a two-stage launch vehicle. As the 
initial point we will take the end of the powered section of the 
first stage. We will assume that the program for this section is 
assigned, more exactly, that is depends on one parameters — the 
maximum angle of attack o m . By selecting a m the direction of the 
tangent to the trajectory at the beginning of the flight of the 
second stage is determined, that is, the angle 0p. 

Depending on the values of the angle of the velocity vector 
9 extra obtained by solving the variational problem In each specific 
case, we will rotate the entire first section of the trajectory so 
that the angle 0p coincides with Op extra and then we will repeat 
the calculation with the revised values of vp, xp, and yp. It 
can be assumed that this method gives a better approximation to 
the extremal program for the entire powered section as a whole. 

By substituting 6u-j_ from Eq. (1.119) into the expression /62 

(1.114) and considering that the variations 6uf and 6yp at the end 
of the powered section are equal to zero, and 6xp and 6yp are 
arbitrary, we get 


pep 




Since wp = vp sin 6^ and up = v^ cos 6p, 


‘g e i=-7 2L =-tg'Pi; 

*11 

— 1; ) 


^sf == 



( 1 . 120 ) 


( 1 . 121 ) 


On the other hand, from system (1.118) it follows that when t * 0, 
we have 


X 

X 


u 

21 




( 1 . 122 ) 
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where tt=-fc tf 


X 1 f=C 1 cos C a sin vtf\ ]j 

A 2 f=C 8 sh |/2 v^ + C 4 ch K 2v/ f ; 

Ag C x v sin vtf f— C 2 v cos vtff; 

X 4 — C 3 |/^v ch y"2v/f — C 4 V 2v sh ]/2v/f . j 


( 1 . 123 ) 


Prom the conditions (1.121) and (1.23) we obtain two equations 
for the four unknown coefficients Ci, C 2 , 03 , and C 4 : 


C,v sin vtff— C 2 v cos vtff— — ; * 

C a y 2v ch'|/2 vtff -f ^4 sh V 3 vtff « ^ 


(1.123' ) 


or, by solving (1.123'), for C 2 and C 3 , we have 


■ C,~ 


+ CjV sin vtf.f 
V cos v*f 

1 — C 4 y 2 v sh \ 2 vtf 
y ~lx ch y 2v(f 


We will find the two equations we lack by using the system 
of equations of motion (1.110). From the first two we get 



'f 

J p sin ydt — S<yz — Wf + Ax* 



where 


Cl : 124 )) 



are the corrections j that 
allow for the nonparal- 
lelity of the field of (1.125) 
gravity and the change 
in acceleration with 
altitude ; 
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C 3 sh y 2 xt + ch y r 2xt 
C 1 cos xt ■+■ Cq sin xt 


( 1 . 126 ) 


9=arctg 

M 


2 = arctg 


The quantities wf and Uf in (1.124) are assigned. 


Considering that Aw and Au are small and act as corrections, 
to the first approximation we can neglect them. Then (1.124) can 
be considered as one of the two equations we lacked for determining 
the arbitrary constants Cq, C 2 , C 3 , and C 4 (or instead of C4, we 
can take tg 0q, since C4 = Cq tg 0q) . The fourth and last equation 
is condition ( 1 . 118 ), in which ul anddwq are the numerator and 
denominator of (1.124). Thus, to the first approximation in 
solving the problem of selecting the program 4 >(t) that gives the 
greatest altitude for the powered section, with assigned values of 
the projections of velocity uf and wf, we must determine the quan- 
tities Cq, C 2 > Co and tg 0q,,by using Eqs. (1.123), (1.124) and 
(1.118). To do this, we can assign several values of Cq (in terms 
of which all the remaining unknowns C2 , C 3 , and C 4 are expressed) 
and for each of them find from (1.124) by selection the 6 q cor- 
responding to it. By then substituting the pairs of Cq and 0q 
we have found into Eq. (1.118), from the series of resulting vi 
we select the value corresponding to the assigned value. Next, we 
c ompute the approximate values of the coordinates x and y based 
on the formula 


(5 ) 

x^x x J rii\t + p/ii j' p cos 9;/$; | | 

6 6 . | 

- i 7 f ! 

= p/Tl j‘(jc 7 sin 9 — | 1 (1.127) 


Let us use the resulting values x(t) and y(t) for computing 
the corrections Au and Aw with which the entire problem is now 
solved in the second approximation. Simultaneously the assigned 
Vq, xq, and yq are revised, as stated above. In addition, by using 
and y-£>, we can, if desired, obtain more exactly the assigned 
direction of the velocity vector at the end of the power section 
(for assigned 0 f = 0 ), and revise wp using the formula 


Wf— 


- y f 


R -I- 


Uf. 


( 1 . 128 ) 


Two approximations are sufficient for the solution to the /64 

problem. The final values of the coordinates are calculated by 
the formulas 


6 l 



Cl. 129) 


JC = x; + a[t + 1 dr\ \ p cos <fUi + AX; 

* V 2 } 

y=y‘ l -lr'w[t+ I 

0 0 *" I 

Here x’^, y T i> u T ]_, and w *l are the following values corrected 
through the first (or second, if required) solution to the problem: 



(1*130) 


To verify these working results, it is useful to integrate 
the equations of motion (1.110) with the above-obtained values of 
4>pr ( t ) . 

We can considerablyr.simplify the problem in practice without 
reducing accuracy if the program <bp r (t) is found for a plane Earth 
and a parallel field of gravitational force, and if in the 
equations of motion (1.110), we take into account all the addi- 
tional terms, regarding them as corrections — certain functions 
of time. 

For the plane Earth, instead of (1.24), we get 


^ p sinarctg -i~ — w f+ 


Sw 


— = tg6i — 




j pcosarctg 4* tg B[ j di — 


Au 


(1.13D 


where C 2 is selected so that Eq. ( 11118 ) is satisfied and Au and 
Aw are determined by (1.125) a as in the first variant of the 
problem. 

If it is assumed that the accleration of the launch vehicle 
n OSO 

p = where n 0 is the initial g-load, and 3 is the relative 
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. tf 

fuel consumption (3 = G/G Q ), the integrals / p sin <}>dt. and 
tf 0 

/ ip cos tj)dt can he written as 


'f 


\p sinarctg (J - + tg ej <lt - -2&&- [-^ In 1 2 (a X)'* + 


~j~ 2ux -j~ b 


b 2 

-f -g— In I I j£- + ft 

' ,1/2 jr > x 1 


'f 


Jo 


| p cosarctg (—■ 4 - tg 0 dt — — n °' g|)C — In 


be W 


2 (cX)W | 2c j ^ 
j: ,v 


'f 


(1.32) 


Here 


X ~<2a; 2 4- £a'4~ c; 

x=a-f^; 

a — 1 /& 2 ; b = — 2ajb\ c = 1 + «‘V^ 2 - 


Based on these formulas, calculations were made for one of the 
possible variants of the two-stage launch vehicle for which the 
dependence of the elements of the powered trajectory section at the 
moment of staging on the angle of the velocity vector, in turn 
dependent on the angle of attack a m , was determined in advance. 

Pigs. 1.7-1. 9 give the altitudes at the end of the powered 
section (the orbital perigee h^) as a function of velocity vp for 
its horizontal direction for three values of the relative satellite 
weight G sa /G (3 = 0.027, O.O 36 , and 0.045 (Go is launch weight). In 
the same three figures are presented the curves of the apogee 
altitude h a (the point most distant from the Earth) and the approxi- 
mate number of satellite revolutions around the Earth during its 
lifetime . 

In the calculation of the insertion section it was assumed that 
3 = 0.005 and no = 1.5 -for the second stage. The reduction in the 
payload in the corresponding variants was compensated by a pro- 
portional increase in the fuel weight, and, therefore, in the 
f lightttime . 

The curves were plotted in terms of the magnitude of velocity 
at the end of the powered trajectory section v^, which the rocket 
attained . 

Analysis of the graphs leads to the conclusion that there is 
an optimal number of satellite orbits (that is, its lifetime) which 
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corresponds to wholly 
specific values and 
h a in each case. For 
the variants under con- 
sideration, the optimal 
orbits are quite elon- 
gated ellipses, and the 
oblongness increases 
with decreasing satellite 
weight . 
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Fig. 1.7. Dependence of h , h a , 
n on vj* (G sa /GQ = 0*027). 
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These calculations 
of the trajectories 
of insertion into orbit 
where 3 « const showed 
that finding the optimal 
insertion program is a laborious 
process. A particularly large 
amount of time was expended in 
computing by the trial and error 
method the integrals appearing 
in Eqs. (1.131) and obtaining the 
required corrections. On the 
other hand, selection of the pro- 
gram and calculations of the in- 
sertion trajectory can be done 
sufficiently accurately, since 
here the optimal characteristics 


Fig. 1.8. Dependence of h^, 
h a , and n on Vf (G^o/Gq - 
= 0 . 036 ). 


7680 y^.m/sec 0 f the orbit are determined or, 
if the orbit is assigned, the 
maximum satellite weight. 
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•r 

j p sin a rc tg (fit + t 0 ) dt — gotf— wf+ A w 
j pcosarctgt^ + to)^— «f-}-Att 


(1*133) 
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To accelerate the calculation 
of the optimal insertion programs, 
we will assume the initial equa- 
tions to be in the dimensionless 
form. 

The initial system of equations 
and the formulas for calculating the 
optimal program (two-stage launch 
vehicle) is as follows: 
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2 . 


*S=® 2 o+«i>; 


a®=2v 2 l t/rfh 


ku= — v a ( xdh 


JC £ ==Jf 0 + i * 0 ^f4- J j/?cosarctg(*ii + t 0 )rfn; 

'f 5 *2 

y^o+'^Vf-f f fife f />sinarctg(6tl+T 0 )rfti— g 0 — ±; 


r 


• 'f 5 , f s 

• ,: f=^+B;%+ j ds| ^cosarctg(6T)+T 0 )dr]-f A *; , A » = 2v 2 j </5 JydH; 

W 2 I * f £ 

jPsinarctg(*t]-j-r 0 )rfTi — g b J^L + A y ; | A*= — v 2 dl^xd^ 


t X f \ 

®f “f « + (for the ease of insertion when 

f 0f = 0); 

The following notation was adopted in these equations: 

P_ z n offo is the launch vehicle' acceleration; 

is the initial g-load 
p! is the fuel consumption coefficient 
b, To : are parameters determining the program 

x 0 =tg: 0 o is the slope of the initial velocity vector with 
respect to the Ox axis 

v 2 _. Jfn . 

R ' ' 

Aa>, A are the additional terms that allow for the effect 
l on velocity of the nonparallelity of the gravi- 
tational force field and the variation in g with 
altitude 


5. 

6 . 
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x- y ■ are the coordinates of the launch vehicle calculated 
without reference to corrections x and y that al- 
^low for the effect of the same parameters as for 
w and y 

x 0 j YOi w 0 j u 0 are the elements of the beginning of the trajec- 
tory section for which the optimal program is 
* determined; the same quantities with the stroke 
refer to the field of revision 

Xf, Yfa Wf, Uf are the elements of the end of the powered section, 

and 

is the altitude of the end of the power section. 


The process of calculation reduces to the solution of Eqs . 

( 1 . 133 ) and ( 1 . 3*0 for the unknown b and tq, with given xo, yQ, 
vq, w^, and u^. 

In the first approximation the equations are solved without 
reference to the correction Aw and Au. Based on the values of uq , 
wg , b, and tq found, let us find the approximate values* off the 
coordinates x and y on the basis of which the corrections Aw, Au, 

Au, and Ax are calculated. In addition, from previously prepared 
graphs for the functions xq, yo, and vq on 0q, characteristic of 
this launch vehicle, let us compute the new revised quantities 
xq 1 a yo 1 3 and vq ' based on the 0Q(t) value found. With these new 
values of the elements of the initial point the calculation is 
repeated, that is, we again select b, tq, uq, and wo- We note 
that if the initially assigned values prove to be unsuccessful,.' it 
is best (If we intend to obtain the horizontal direction -of "the /69 

velocity vector Vf) to revise w f by means of Eq. (1.128). The 
process of approximation can be repeated several times, depending 
on the required accuracy. 

Finally, after carrying out all the revisions the altitude 
of the powered section h^ is computed. 

Let us bring the integrals appearing in Eq. (1.133) to the 
dimensionless form. To do this, we introduce the dimensionless 
time $ = t/tf. 

The integrals become: 

f p sinarctg(A/-j-T 0 )^ = « 0 g- 0 ^j, 

6 

sinarctg(6'$ + t 0 ) ^ 

1 — 
o 


(1.35) 
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? 


p cosarctg {bt -f t 0 ) dt=n (i g ( ^i 2 , 


r _ e cosarctg^'fc + xp) ^ 
1— a£ 


j ^ \psmarctg(bt+T ( ,)dt=n 6 g lt t‘/ > , 


X 71 

A= ^ -gj nar ctg(*^4-To) 

6 o 1 — 




'-f ? 

j ti j p cosarctg {bt -f t 0 ) dt=n 0 gJ*f A , 




M 

0 0 


cosa rctg (b'j + t 0 ) 
1 — a£ 


d\\ 


[y^=--« 0 4 r 4 3 / 61 


ooo 


sinarctg(6'$ 4- Tp) 
l~a£ 


dV> 


*£ 

j xdt=n 0 g^ 9 , 


/,= \ rfE ( rfti C ;°?»L^tf: . e±M rf6 ; 

OOO 


1 — flt 


r*f 


j ydt^n^y ^ 


ff o 


i i ^ 


/,= J <*« r sina f ctg(^ g + T 0 ) J; . 

n f\ a it ^ ^5 


0 o 6 o 


'f* 6 
j d\ j xdt—n^g^l 
6 o 


8 > 




I I n (i 

a 6 o i 


I I 


= J d\ | (in jd|i j V cosarctg («-£ + !„) ^ 


(1.36) 


(1.37) 


( 1 . 38 ) 


(1.39) 


/70 


(1.40) 


(1 . 4l) 


(1.42) 
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Here: a is the fuel filling factor (a = G m k/^01> where Gq-, and 

ni]_ is the initial weight and the weight of fuel in 
the second stage), and 

b* = btf (the stroke is omitted in the following). 

Let us rewrite Eqs. (1.133) with reference to the expressions 
obtained for the integrals !]_, ...» Iq: 


1 . 


2 . 


3. 


4. 


• _ ffi, o _ 

irt " 


“0 


— n n g Q tfI i 4- f— Aw 

— 2 + — Au 

®5=®o+«o; 


A® — 2v 2 f ydt ~ 2v 2 go ^ b 


'f 


AW= — v 2 ^ — + "ogo^e i 


Ay — 2v a 


•*f— + «(/f + w ogo^4 

/ 2 

?f= y o 4- «Vf4- nogo'ps - go -f-; 

x 4 Q -\- u'^-\r n^g Q ^{ A -\- A*; 

tjf~ y' 0 4 »p!?o^ 8 — g'o + Ay ; 

/ 2 / 3 1 
■*o — r u o 4" n ogo^ 8 J 5 

-t^- 4- ™o~“ — -—--Wiogo^j; 


A* 


y 0 
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5. 


6 . 


Wf 5 = ~tt£- 


//f-yff 


/? 4 - 
2 /? ‘ 


*f 


The initial point (xg, yg) must be selected beyond the limits 
of the dense atmospheric layers, since as indicated above, the 
program of the powered section is selected in the dense layers 
independently of the solution of the variational problem. 

This method of calculation can be used not only for two-stage, 
but even three-stage rockets. By a stage we mean sections -of a 
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